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2 : ; Um umfangreiche lineare Gleichungssysteme auflösen zu können, wird die gegebene Koeffizienienmatrix 
Br = in Untermatrizen aufgeteilt. Mit dem unterteilien modernisierien Ga ußschen Algorithmus werden die 
Be 3; Unbekannten gruppenweise eliminiert. Durch Binschalien von Zwischenkontrollen wird die Rechensicherheit 
ee ‚erhöht. Es ist die Möglichkeit gegeben, die Arbeit auf mehrere Rechenkräfte gleichzeitig aufzuteilen. Die 


N 
\ 


Vorgehensweise läßt sich auf die Berechnung der Kehrmatrix ausdehnen. 


To solve extensive systems of linear equations, the coefficient matrix: is divided into sub-matrices. Then, 
the unknowns are eliminated, in groups, by aid of the sub-divided abridged Gauss algorithm. Thus’'it is 
possible to accomplish the solution by several reckoning. assisianis simultaneously. The security of reckoning 
is increased, if intermediate controls are inserted. The method is applicable, 100, to the-computation of the 
inverse matrix. h : 


Kacı 


Pour V’Evaluation de syst&mes d’&quations lineaires tres &iendus la matrice donnee de coefficienis est 
diviste en sousmalrices. Les inconnues sont &limintes par groupes a l’aide de l’algorithme divise et modernise 
de Gauss. La süret& de calcul esi augmenite par des contröles intermediaires. Il est possible de partager 
le travail calculatoire entre plusieurs calculateurs simulianement. Le procede peut &tre &tendu pour V’&valu- 
ation de la matrice inverse. 

Ara pemenns 601BmHX CHcTeM AuHeüHsx ypaBueuni anna Marpuna koabhunmenTtoR 
pasönsaerca Ha HOAMATPHIIEI. Heussecrusie HCKAIIOYAPTCH HO TPYTDAM IPH HOMOINH HOAPA3- 
= AeleHHOTO YCOBePINIEHCTBOBAHHOTO anropndma T’aycca. Hane:kHoCTs pacyeTa MOBEITAETCH 
= HPOMe:KJTOYHEIMH KOHTponAaMmnm. Pabora Mo:KeT ÖbBITB pacıperereHa ONHOBPeMEHHO H& He- 
= CKONBKO Bhlyncauresneä. Meron aonyckaer pacıpocrpaHeHunue Ha BBIyHcHeHue O6bpaTHoh 

.. M&TPuIBI. 

3 Einleitung 
- Bei der Auflösung linearer Gleichungssysteme mit einer großen Anzahl von Unbekannten 
_ sind zur rationellen Durchführung der Zahlenrechnung oft andere Gesichtspunkte zu beachten 
_ als bei Gleichungen mit nur wenigen Unbekannten. In der vorliegenden Arbeit wird eine Vor- 
_ gehensweise beschrieben, die den besonderen Anforderungen umfangreicher Systeme angepaßt 
_ ist und die sich praktisch bereits gut bewährt hat. Zur Beschreibung bedienen wir uns der 
_ vereinfachenden Matrizenschreibweise. | er, 
Dem Verfahren liegt der im folgenden als bekannt vorausgesetzte sog. modernisierte 
_ Gaußsche Algorithmus zugrunde (weiterhin abgekürzt mit m. G.A. bezeichnet), eine durch 
zweckmäßige Anordnung und Ausnutzung der Rechenmaschine herbeigeführte Vereinfachung‘ 
es gewöhnlichen Algorithmus. Die Methode ist ausführlich beschrieben bei A.Walther und 
H. Unger [16], für Normalgleichungen (symmetrisch positiv definite Matrix) bei E.Andersen [2] 
‘von dort ist auch die Bezeichnung m.G.A. übernommen), ferner bei R.Zurmühl [21], [22], 
. 248263, wo besonders auf den Vorteil des Einsatzes der Rechenmaschine hingewiesen wird, 
as wohl zuerst bei Cholesky [4] und auch schon bei Andersen [2] geschieht. Mit der gleichen 
frage befaßt sich Th. Banachiewicz in einer Reihe von Arbeiten [3] (vgl. auch [9], S.16). Der 
A. findet sich’in der älteren Literatur auch als ‚„abgekürzter Gaußscher Algorithmus‘‘, 
h ohne Angabe der durch die Rechenmaschine zu erzielenden Vereinfachungen (z.B. [10], 
ildung der Endgleichungen ohne Zwischenglieder). Re 
sim m.G.A. sind Ausdrücke von der Form 
SE ab+cd+tef+-- 
, wobei man die Teilprodukte in der Rechenmaschine selbsttätig zum Gesamtausdruck 
äßt. Dieser Vorgang kann nun bei Gleichungen mit sehr vielen Unbekannten praktisch. 
liebig viele Teilprodukte ausgedehnt werden, ohne daß die Sicherheit der Rechnung 
beeinträchtigt wird. Diese Überlegung führt dazu, die Elimination der Unbekannten 
Schritte aufzuteilen und die Rechensicherheit durch Einschaltung von Zwischen- 
‘erhöhen. Hierdurch tritt keine nennenswerte Vermehrung der Rechengänge ein. 
er is es durch eine solche Unterteilung möglich, die Arbeit auf mehrere gleichzeitig 
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anzusetzende Rechenkräfte aufzuteilen und’ dadurch eine ganz wesentliche Herabsetzung der 
Arbeitszeit herbeizuführen. gr i 
Hierzu wird die gegebene Matrix in eine Anzahl von Untermatrizen aufgeteilt. Man arbeitet 
in Teilschritten und eliminiert jeweils mehrere Unbekannte gleichzeitig. Das Vorgehen entspricht 
dem gewöhnlichen Gaußschen Algorithmus, wenn man als Elemente Untermatrizen anstatt 
gewöhnlicher Zahlen nimmt. Dabei tritt an die Stelle der Divisionen die Bildung von inversen 
Matrizen. Es zeigt sich aber, daß die Bildung der Kehrmatrizen zu den Untermatrizen, die 


einen Umweg darstellen würde, nicht explizit notwendig ist. Die Vorgehensweise bezeichnen wir: 


als unterteilten m.G. A. | 

Das gleichzeitige Entfernen mehrerer Unbekannter istfrüher schon von R.Mehmke[l2], [13] 
behandelc worden. Wir zeigen hier, daß sich diese Aufgabe zweckmäßig mit dem unterteilten 
m.G.A. erledigen läßt. : 

Im folgenden wird keine spezielle Gleichungsmatrix zugrunde gelegt. Vielmehr soll jede 
beliebige quadratische Matrix mit von null verschiedener Determinante auch bei sehr großer 
Reihenzahl bearbeitet werden können. Die Elemente der Matrix werden hier als reell ange- 
nommen, doch läßt sich alles ohne weiteres auch auf den Fall komplexer Gleichungssysteme 
übertragen. Bezüglich einer Umschreibung eines komplexen auf ein reelles Gleichungssystem 
mit 2» Unbekannten vgl. etwa R. Zurmühl [22], S. 269. 


1. &außscher Algorithmus mit Untermatrizen 
Gegeben sei das inhomogene lineare Gleichungssystem 
Yı=e. "mit Ve), vIeba san BR 


mit nichtverschwindender Koeffizientendeterminante det Y++0. Gesucht ist der Vektor x für 
eine gegebene rechte Seite a. Zur Behandlung umfangreicher Systeme werde (1) folgendermaßen 
aufgespalten ?): 


r 


» 
r=1 
Dabei sind die U quadratische Untermatrizen mit einer Zeilenzahl 07, wobei 


+%+'+%=n. 


DT en 2 0)4 


Die übrige Aufteilung der Matrix X ergibt sich dadurch zwangsläufig; die Matrizen X) haben 


co, Zeilen und o, Spalten und seien kurz als (o,, o,)-Matrizen bezeichnet. ‘Wir setzen für das 
folgende die in der Diagonale im Verlaufe der Rechnung erscheinenden quadratischen Matrizen AY 


als nicht singulär voraus, was nötigenfalls durch Umstellung der Gleichungen oder Unbekannten 
erreicht werden kann. 


Auf (2) kann man formal den Gaußschen Algorithmus anwenden, wobei anstatt mit 
Zahlen mit Untermatrizen zu operieren ist. Man multipliziert demgemäß die erste der Gl.(2) 


a a1 : i n 
‚mit Wi und dann mit — A} vor und addiert die so abgeänderte erste Gleichung zu der 
ten, 2=2,3,..., p. Dadurch werden die o, ersten Unbekannten eliminiert und man erhält 


g ar®) (8) \ 
Zr u=0 mit Lan seen DT EEE TRER (3). 


Die Untermatrizen un entstehen dabei durch die Vorschr 


iften des Eliminati S 
wandt auf Untermatrizen a = 


(e) ! 6) leer 
Kr AN ULRU nn ee Ehe 


Auf Grund der getroffenen Unterteilung sind die X auch (07, o,) Matrizen. . 
Aus dem System (3) werden durch Elimination (Bildung von un 


entfernt. Diese Vorgehensweise wird fortgesetzt, bis man das Gleichungssystem: 


erhält, in dem nur noch o, Unbekannte 


dann die Aufrechnung aeg ganzen Systems erfolgen, indem rückwärts eingesetzt wird. Mit Hilfe 
der Kehrmatrizen (rn) können die unbekannte | 
P—23 ..., 1 bestimmt werden. 


2) Der obere Index in A) und al) 


) weitere o, Unbekannte 


U LT ee (5)52 
auftreten. Nach Lösung des Gleichungssystems (5) kann R 


n einspaltigen Teilmatrizen x, fürr=p—1, 


n E ; soll zur Unterscheidun S ’ A 
daß es sich um die gegebenen Ausgangswerte handelt. lung’ von späteren Matrizen darauf en 
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2. Zweekmäßige Durchführung des Algorithmus mit Untermatrizen 
Der Nachteil der soeben geschilderten Vorgehensweise besteht darin, daß bei jeder Teil- 


elimination von o, Unbekannten zunächst die Kehrmatrix A) " und dann gemäß (4) das 
Produkt dreier Matrizen gebildet werden muß, was einer wesentlichen Erhöhung der Rechen- 
gänge gegenüber dem gewöhnlichen Algorithmus gleichkommt. 

Dies läßt sich vermeiden, wenn man auch auf die Teilelimination den m.G.A. anwendet 
und jedesmal nur o, Unbekannte eliminiert. Hierdurch wird natürlich die Schreibarbeit gegen- 
über dem m.G.A. in einem Schritt erhöht. Dies aber wird durch erhöhte Sicherheit, durch 
Anwendung von Zwischenkontrollen und vor allem durch die Möglichkeit der Arbeitsaufteilung 
auf mehrere Rechenkräfte wieder wettgemacht. 


Der m.G.A. bedeutet im Matrizenkalkül die Aufspaltung einer quadratischen Matrix in 


das Produkt zweier Dreiecksmatrizen (vgl. [16], [2], S.11; [22], S.250). Angewandt auf X 
heißt dies: 


(1) 
U DR a ARTE (6), 
ausführlich: = 
1 1 
NE as, ET FREE Ui Bene ae 61,0, 
1,0)" @ (a) De DEN SUR Oi Ca 

IN EIER 20, Iaı 22 ‚or 
eo... x SS te 63) 
Ga,,1 A022 "Her. (ig; ‚0, Qa,1 Qa,2 rer. 1) ER Co,,0; 


QD,, ist eine untere (N) und @,, eine obere (N) Dreiecksmatrix. Die Rechenvorschrift zur Er- 
mittlung der Elemente in Q,, und @,, ist gerade der m.G.A. 


Aus (6) folgt für die Kehrmatrix von An: 


= 2 2 i 
a Ce (7). 
Zur Durchführung der Vorschrift (4) benötigt man also die Matrizen 
| Ar Ca = Aıı und ee) 
womit (4) sich vereinfacht zu 
an = an EN N ee a ner er (4a). 


Schema zur Auflösung umfangreicher Gleichungssysteme 
(Aufteilung in neun quadratische Untermatrizen mit je o Zeilen, 30=n.) 


3 
Gleichungssystem: 2% An > a mit D=—1,72, 3, 


nt 


Sirenen] 
elgelyeil 


EEE Bee] 


Bild. 
1 * 
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i ot) 
Die Untermatrizen Q,ı und €, erhält man aber einfach durch Anwendung des m. G.A. auf aa 


bzw. A) mit den bereits bekannten Dreiecksmatrizen Yyı und &,, gemäß: 
a) 


r x 
Dr &,, = Ui und Du Cr Nr Per ee 2.788); 
Zur Erläuterung sei ein Schema (Bild 1)gebracht, in dem in 3.3 9 quadratische A 
gleicher Reihenzahl o aufgeteilt sei. Bei’der ersten Elimination werden o,—=o Unbekannte en i 
fernt. Die Durchführung dieses Prozesses mit den UNISRE EI Dun Dzı D,ı und E,,, Ca, Cs 
und das Restsystem mit den Untermatrizen A, Xi , Ass und X5s ist im Schema angegeben, 
während die aufgeteilte ursprüngliche Matrix nicht aufgeführt ist. = 
Die neuen rechten Seiten werden in den Matrizen Q,;, 2, AR als zusätzliche (o+1)-te 
Spalte mitgeführt. 


Elimination der Unbekannten %, bis u. 
F ; di), 
Berechnung der Elemente in &,,, Cs, &; nach A He: 
Gar — Fir — hear —'— mir für I=1l,2..,0; (D). 
k=i,1+1,..,,0+1..,20,204+1,..,0, n+1 
Berechnung der Elemente in Qu; Day, Dı nach Dyı HM: 
Gk = [ — qui De ET ee nr 1,21 %-1,8]: Gr. - ürsk= 1 2 (II). 
: i=k+1,..,0,0+1.,,20,20-+1..53%2; 1 
A Zuerst Bestimmung der Matrizen Q,, und @&,, wie folgt: 
Be Auswertung von (I) für «=1, k=1 bis o, dann (II) für k=1, i=1 bis o, dann wieder (I) 
: 7 - fürs=2, k=2bis o, (II) für k=2, «=2 bis o usw.°) Nach Ermittlung von Q,, und @,, können 
ge er dann Q,,, DO und GC, C,, unabhängig voneinander bestimmt werden. ; 


E ei | Berechnung der Elemente in AP für ,r—=23 nach UN — Dr rn Eır: 
og } 

@ = ar — qı C18 7 1802. IR Nieler: SIE 
Mn i=0+1,0+23,..,n; k=0o+1,0+2%,..,nn-+1 

EB; Kontrollen: 


s und ?sind die Spalten, bzw. Zeilen für die Kontrolle durch Zeilensummen, bzw. Spalten- 
'  summen. 


Die Elimination der Unbekannten x,;ı bis x. erfolgt beim nächsten Schritt. Schließlich 
wird AS} in zwei Dreiecksmatrizen aufgespalten. - . 
Die Aufrechnung der Unbekannten erfolgt nach: 


6: Ci,n+1 Ci,n %n — in In-1 090 Girl Bir i=nn—l..„1 (DV). j , - 
Aus den Herleitungen ergibt sich: ae PER Dr r 
dtX=|U=|E&ıl’|&sl | -»--- :- N 


Damit kann auch die Auswertung umfangreicher-Determinanten auf diese Weise vorgenommen 
werden. Bereits A.Walther weist darauf hin [18], daß die Auswertung von Determinanten 
' ‚nach der Regel von Chiö [19] auf den gewöhnlichen Algorithmus hinausläuft. Die Aufspaltung 
der Matrix X in zwei Dreiecksmatrizen O und € führt dann auf die Verwendung des m. G.A. und 
für umfangreiche Determinanten auf den unterteilten m.G.A. 

» k 


3. Praktische Durchführung 2 a 
: ‚ Für umfangreiche Rechnungen sind laufende Kontrollen unerläßlich. Die beschriebene 
Vorgehensweise ergibt die Möglichkeit, jede der auftretenden Untermatrizen auf ihre Richtigkeit 
zu. prüfen (notwendige Bedingung). Im Schema sind dazu Kontrollen durch Spalten- und 

-  Zeilensummen vorgesehen, die wahlweise oder zusammen benutzt werden können. ? 


r _ Bei der Rechnung ist es unter Umständen zweckmäßig, bei der Aufspaltung von x 
‚spalten- oder zeilenweise vorzugehen. Man kann dann im Falle eines betragsmäßig kleinen Haupt- 
diagonalgliedes einfach eine Umstellung vornehmen. Hat man schon ko Unbekannte entfernt, 


dann kann man die Umstellung in der restlichen Rechnung allein vornehmen, ohne die vorherige er 
% een ö j 1 aa Pr ? Er Be as 
ar.  .®) Vgl. hierzu auch Abschnitt 3, a SH er: gr 


Tg 
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. Rechnung abändern zu müssen. Lediglich beim Rückwärtseinsetzen hat man auf die Reihen- 
folge der“Unbekannten zu achten. 
Neben erhöhter Übersichtlichkeit und Kontrollmöglichkeit ergibt sich bei der Unter- 
matrizenmethode als wesentlicher Vorteil die Arbeitsaufteilung auf mehrere gleichzeitig einzu- 
_ setzende Rechenkräfte, nachdem erst eine gewisse Anfangsrechnung, nämlich die Berechnung der 
beiden Dreiecksmatrizen Q,, und @,, durchgeführt worden ist. Hierdurch-kann die für die 
Rechnung benötigte Gesamtzeit wesentlich herabgedrückt werden, wobei selbstverständlich die 
Gesamtzahl der Operationen die gleiche bleibt. Die zweckentsprechende Art der Aufteilung-wird 
von Fall zu Fall leicht gefunden werden können, so daß sich nähere Ausführungen erübrigen. 


4. Anzahl der Operationen und Zeitbedarf 


Hier sei wieder der Einfachheit halber angenommen, die Matrix sei aufteilbar in p2 quadra- 
tische Untermatrizen gleicher Zeilenzahl o. Dann sind, wenn wir auch die Bildung von 9, © 
als Berechnung einer Matrix auffassen, insgesamt 


- 1 
E S=srP+NDeP+1) Be har a re era ee (10) 
E Untermatrizen zu berechnen. Davon erfordern?) (ohne Kontrollen und ohne rechte Seite) 
- | 3 2 
N ee » Untermatrizen 5 5 + 7 Multiplikationen, 
a 2 | &e ERS Ku ja 
a p p 9 3) . ’ ’ 
1 
5 PP—-U@P—1) » o® “ 
i Insgesamt erfordert die Aufspaltung von X in zwei Dreiecksmatrizen: 
u 
e- : BOT v0° om nn ae Son 
7 | —£@2._P% 1 P2 2” ._9 1% Multiplikationen . 2... 11 
E Z, 3 5 + EL PT -4 5 Multiplikationen (11) 
4 mit po=n, wie beim gewöhnlichen Gaußschen Algorithmus. 
= Von den $ Untermatrizen brauchen aber nur (im Falle oe >1, p<n) 


EEE 
\ 


TERTINU REN 


re a (12) 
nacheinander berechnet zu werden. Hiervon erfordern 


MEEREON 0: KG SERn en 
p Untermatrizen Sn a, 6 Multiplikationen, 


1 erde | 
p == „ 92 5) „ 3 
p nt an 0? „ D 
also insgesamt: ; > | 
£ Ga = Se 3) —5 Op z Multiplikationen . . 2... (13). 


Das Verhältnis a _ q gibt näherungsweise einen Anhaltspunkt für den Bruchteil der 
Zeit, auf den man durch Arbeitsaufteilung herunterkommen kann und damit auch für die Zahl 
der Rechenkräfte. A ES 

Dee Es kann unter Umständen zweckmäßig sein, als Matrix Wii eine Matrix mit wenig Spalten 

_ zu wählen, um möglichst schnell mehrere Rechenkräfte einsetzen zu können. 

Er Um eine gleichmäßige Aufteilung zu erzielen, kann es notwendig werden, einzelne der 
- Untermatrizen noch weiter aufzuteilen. Fe: ' EB 
Eine Aufteilung in vier Untermatrizen bringt naturgemäß nur wenig Zeitgewinn, da mit 

= j o%-. 0? N . 
| Ausnahme der beiden Matrizen D,,, Cıs (@-3) Muliplikationen 2 Ber x 

rechnet werden muß. Ebenso kann die Arbeitszeit bei Matrizen, die nur aul und in der Nähe 

x Hauptdiagonale besetzt sind, oft nicht verringert werden; immer bleibt aber der Vorteil der 

ersichtlichen Anordnung und der guten Kontrollmöglichkeit erhalten. 


4 erden air die Multiplikationen als wesentliche Operationszahl betrachtet. 
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Beim Rückwärtseinsetzen können entsprechende Unterteilungen vorgenommen werden, 
nachdem zunächst Wat» — ap, gelöst worden ist. 


5. Berechnung der Kehrmatrix 


Die Aufstellung der Kehrmatrix X-1 zur gegebenen Matrix X ist bekanntlich sinnvoll, wenn 
das Gleichungssystem (1) nicht für eine gegebene Seite, sondern für mehrere, u.U. noch frei- 
bleibende rechte Seiten zu lösen ist. ? 

Zunächst kann man fragen, bei welcher Anzahl k rechter Seiten sich überhaupt die Be- 
rechnung der Kehrmatrix lohnt. Dazu ist festzustellen: Unabhängig von k ist die Berechnung 
der Kehrmatrix immer ein Umweg, da sie größenordnungsmäßig 2 n?/3 Multiplikationen mehr 
erfordert als die direkte Rechnung. Denn es müssen ja nicht nur die Operationen zur Bildung 
von A-! gezählt werden, sondern auch diejenigen, die notwendig sind, um X" mit den rechten 
Seiten zu multiplizieren. Die zusätzlichen 2 n°/3 Multiplikationen sind gerade diejenigen, die 
notwendig sind, um die Kehrmatrix zu berechnen, nachdem die Aufspaltung.von X in DE 
durchgeführt ist (vgl. [16]). 

Damit soll natürlich nicht in Abrede gestellt werden, daß in gewissen Fällen die Berechnung 
der Kehrmatrix Vorteile bringt. Insbesondere ist der Rechnungsgang für eine neue rechte Seite 
bei Kenntnis der Kehrmatrix übersichtlicher. 

Um das in 2. beschriebene Verfahren zur Berechnung der Kehrmatrix A! anwenden zu 
können, gehen wir wie üblich von dem Gleichungssystem 


DR x m E mit x —— An ee cl RE ER (14) 


aus. Das Schema wird mit den n Spalten der Einheitsmatrix € als rechter Seite durchgerechnet, 
wobei die Aufteilung in Untermatrizen sich aus der bei X vorgenommenen ergibt, und wo @,, 
mtr=p+19+23..,.2p für r>(p-+L) Nullmatrizen und @,,»;: untere Dreiecks- 
matrizen (N) darstellen. Die Aufrechnung der Unbekannten für die n rechten Seiten liefert die 
gesuchte Kehrmatrix, wobei die gleiche Unterteilung beibehalten wird. 

Eine andere Vorgehensweise ist die folgende, bei der wieder eine Unterteilung in Unter- 
matrizen vorgenommen werden kann. Als Ergebnis des in 2. beschriebenen Verfahrens erhält 
man — wenn man ohne rechte Seite arbeitet — 


EN ee ee a A (15), 
wobei Q und € Dreiecksmatrizen sind mit 
D= (A) und Ellen i=T ar Per (16), 
und mit 
Dr —=0 juris 27.31 ZuBu Er=derSite 
Durch Auflösung der beiden gestaffelten Gleichungssysteme 


C>E md Du Eee (17), 


erhält man die beiden Kehrmatrizen €! und Q-!, aus denen sich dann durch Multiplikation 
Be — 19-1 ergibt. Wir schreiben die beiden Kehrmatrizen, die wieder releckanies&£ vom 
gleichen Typus wie € bzw. O, sind, in der Form?): 


O=(&,), (G,=0 für I>9), Sm=(&), (&,=0 für I<r)... (18). 
Nur‘die Hauptdiagonalmatrizen Q,, und €, sind wirkliche Kehrmatrizen, { 


Gr Re De . (19). 


Dabei benötigt man zur Berechnung von Q,, die Kenntnis von Q 
gekehrt bei C&,, muß €&,,, &,_1,r°'', Cy41,r bereits bekannt sein. 


Die Unterteilung der Matrizen Q und € und damit auch von O-! und G-1 jst durch die 


ursprüngliche Untermatrizeneinteilung der Matrix festgelest ießli ibt si i 
res g gelegt. Schließlich ergibt sich die Berech- » 


rr>9 DOr+i,r bis Dur und um- 


an De . 
Y-A1— (%,,) u | 2 Cm ön.) En a ae ea (22) Zu 
wobei jede Untermatrix X,, wieder für sich berechnet werden kann. = 


5 Di . . ° , 
En) ie Bezeichnung gibt zu keinerlei Verwechslungen Anlaß, da alle Elemente in Y als reell vorausgesetzt 
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Von Interesse ist noch der Fall p=2, Unterteilung in vier Untermatrizen, Die Kehrmatrix 
ergibt sich in der Form: 


yX-1 Ip, 5 Di ar Ca Dan» PR >) 
Ca Dgı 9 Ca Das 


Schreibt man diesen Ausdruck so um, daß man die vier Untermatrizen in (21) durch u und 
die entsprechenden Kehrmatrizen ausdrückt, so erhält man: 


_1 1 = “ = v 
ae u, tee 
ee, up” Be) 
wobei AS? nach (4) erklärt ist. Diese Beziehung stellt das Verfahren von Boltz dar [6] (vgl. 
auch [9]), das schon von R.Lohan [11] in die Matrizensprache übersetzt worden ist und das auch 
leicht aus einer Arbeit von J. Schur [15] gefolgert werden kann. (Vgl.hierzu auch E.Bodewig [5] 
und Frazer-Duncan-Collar [8].) 
Die Relation (22) erweist sich als ein Spezialfall der allgemeinen Darstellung (20) einer 
Kehrmatrix mit Untermatrizen. Mit den eingeführten Untermatrizen erhält man eine bequemere 


und übersichtlichere Darstellungsform. 
Man erkennt, daß allen Beziehungen der Gaußsche Algorithmus zugrunde liegt und alle 


‘ wesentlichen Zusammenhänge hergeleitet werden können, indem man die Elemente _durch 


Untermatrizen ersetzt. 


6. Symmetrische Matrizen 


Bei symmetrischen Matrizen treten beträchtliche Vereinfachungen auf. Behält man die 
bisherige Methode bei, dann können die Elemente der Untermatrizen Q,, einfach aus den zur 
Hauptdiagonale spiegelbildlich gelegenen Elementen der @,, bestimmt werden. 


dir = Cri : Okk RE IE N ERTERRE ea (23). 


Bei Berechnung der Kehrmatrix brauchen nur die auf und oberhalb der Hauptdiagonale gelegenen 
Elemente berechnet zu werden, da mit X auch X! symmetrisch ist. 

Schließlich kann man auch die Vorgehensweise voll symmetrisieren, d.h. man rechnet nur 
in- und oberhalb der Hauptdiagonale, indem man die Methode von Cholesky [4] (vgl. [9], S.21 
und [22], S.260) heranzieht. Dabei treten in der Hauptdiagonale die Quadratwurzeln der Haupt- 
diagonalelemente auf, durch die jede Zeile dividiert wird. Die Matrizen Q,, werden überflüssig. 


7. Zusammenhang mit weiteren Matrizenoperationen 
Die im Schema von 2. nach Elimination der o ersten Unbekannten nach Gleichung (III) 
erhaltenen Elemente a (i, k> o-+ 1) können als (o-H1)-reihige Superdeterminanten zu [Ar], 
dividiert durch |XS1, dargestellt werden. 


dıı Se ee te ee AUco 41% 
daı Beine she ee A3o 42% 
h She een ole let el d en men nern H P 
Ber . a — N Pl ee a (24), 
(oi EEE Ico Goch 
Al meer en io Air 


‚wie schon in-einer Arbeit von A.Walther [17] näher gezeigt ist. Aus (24) folgt dann, was gleich- 
falls in [17] erwähnt ist, eine einfache Herleitung des Determinantensatzes von Sylvester, den 
wir jetzt in Anlehnung an unser Schema so formulieren können: 


(2) ar(2) 
A 2: 
Asa Ass 


Denn die 2. Determinante links ist ja gerade die Determinante der durch |? | dividierten Super- 


Iaf?]- 


> ® 1 
determinanten zu bien ; 


‚Hier ordnet sich auch die sog. „‚beschleunigte Elimination‘ von R.Mehmke [12], [13] 


unter, bei der ebenfalls mehrere Unbekannte gleichzeitig eliminiert werden. Während jedoch bei 
_ Mehmke hierbei die Kenntnis der algebraischen Komplemente erforderlich war und schon aus 
diesem Grunde praktisch nicht über o=2 hinausgegangen wird, ermöglicht der m.G.A. die 


er r. “ 


AUR EN  1 e 


EM 


# 
j 


- [8] Frazer, R. A., W.J.Duncan u. A.R.Collar: Elementar 
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gleichzeitige Elimination von beliebig vielen Unbekannten mit etwa gleichem Aufwand, jedoch 


ohne explizite Kenntnis der Komplemente. Schon G.Worch [20] bemerkte, daß der Ausdruck 
x hleunigt““ nicht zutreffend ist. Sa 
altken [1] gibt eine Methode zur Berechnung der reziproken Bilinearformen _ 


EV en 


it Hi öhnli i 1 die nichtsinguläre . 
mit Hilfe des gewöhnlichen Gaußschen Algorithmus. Dabei ist gegeben 
Formmatrix Sn die im allgemeinen nicht quadratischen Matrizen ©, und &,. Aus einem 
Vergleich mit (4).ersieht man, daß der Ausdruck (26), wenn man vom Vorzeichen absieht, nach 
der Vorschrift des m.G. A. berechnet werden kann, indem man £ 


1 a) a 
A — 9, Der — G;; Azı >= &,, Aaa ee 0 


setzt. 3 
Auch die Auflösung des linearen Gleichungssystems (1) kann in dieser Form durchgeführt 
werden mit & = — E nach nebenstehendem Schema: Se . Dies läuft im wesentlichen 


auf das Verfahren von Gauß-Jordan [10] hinaus und ist nur in Ausnahmefällen empfehlens- 
wert, da ein Umweg. 


8. Schlußbemerkungen 
. » Die Untermatrizenmethode bietet die Möglichkeit einer Unterteilung derart, daß jeder 
„Block“ in sich fertig gerechnet und kontrolliert werden kann, bevor er weiter benützt wird. 
Die Aufteilung bietet dadurch weiter die Möglichkeit, mehrere Rechenkräfte gleichzeitig ein- 
zusetzen. 

Bei linearen Gleichungssystemen wird die’Zahl der notwendigen Multiplikations-Rechen- 
gänge größenordnungsmäßig durch n/3 bei einer rechten Seite und durch n® bei einer Kehrmatrix 
im wesentlichen festgelegt. Kleine Abweichungen, die keine Rolle spielen, können durch ver- 
schiedenartige Anordnung auftreten. 

Heute können mit gewöhnlichen Rechenmaschinen Gleichungssysteme mit 100 bis 200 Un- 
bekannten ohne weiteres gelöst werden. Weitere Fortschritte sind möglicherweise durch den 
Einsatz moderner Rechenautomaten zu erwarten. 

Es sei noch auf die Möglichkeit einer Verbesserung der Ergebnisse durch Nachiteration 
hingewiesen, falls die erste Rechnung infolge der natürlichen Begrenzung der Stellenzahl unge- 
nügende Genauigkeit liefert. (Vgl. hierzu [16], [22] S. 2832—285 und [14].) 
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| Maximalkorrelation und Korrelationsspektrum 
) Von Hans @ebelein in Landshut (Ndb.) 


Die Theorie der Maximalkorrelation liefert für den Fall einer Korrelationstabelle mit k Zeilen und 
Spalten (k—1) ausgezeichnete Korrelationswerte, die das Korrelationsspektrum, darstellen. Bei einer zwei 
parametrigen, konlinuierlichen Verteilung umfaßt dieses Spektrum sogar im allgemeinen unendlich viele 
Werte, z. B. bei der zweiparametrigen Gaußschen Verteilung die aufsteigenden Potenzen r% des gewöhnlichen 
Korrelationskoeffizienten r. Umgekehrt gibi es aber nicht zu jedem vorgegebenen Spektrum Häufigkeits- 
verteilungen, sondern die Korrelationswerte müssen auf gewisse Bereiche beschränkt sein. 


The theory of maximum correlation provides (k—1) exceptional correlation coefficients in case of a 


: table with k rows and columns, these values constituling the correlation spectrum. T'his correlation spectrum 


Z in general contains an infinite multitude of values if the distribution in question is continuous, e. g. the 
correlation spectrum of the bivariate normal distribution is given by the powers r" of the ordinary correlation 
coefficient r. What concerns the inverse problem, it is not possible to choose a correlalion specirum arbi- 
trarily but its set of values is restricted to certain areas. 


je La theorie de la corr&lation maximale fournit (k—1) valeurs exceptionelles du co£fficient de corr£lation 
Z s’il est la matiere d’un tableau avec k series et colonnes, et ces valeurs forment le spectre de correlation. Ce 
spectre gen£ralement est compose d’une multitude infinie de valeurs si la. distribution est. conlinue, par 
exemple le spectre de correlation de la distribution normale & deus variables contient les puissances r" du 
co£fficient ordinaire r de corr&lation. Pour le probleme inverse il n’est pas possible de choisir arbitrairement 
les valeurs de la spectre de corr£lation, mais ces valeurs doiveni &ire enlour&es dans desr&gions characleristiques. 


Teopnua MakcuMasıbHof KOppenaumn HaeT B CIHyyae KOPPenAnmuoHHON TAOIHIBI C k CTPOKAMH 
a cronönamn (k—1)ompeperteHHoe KOPPpeIANHOHHOE BHAYEHHE, HPEACTABNAINEE KOPPeJIAIHoH- 
Hs cıertp. Ilpu nsymapamerpnueckom, HempepsIBHOM pacnperenennn 3TOT CIIeKTP, B 
oÖINeM CAy4a®, CONepiKUT Narke beckKOHeYHO MHOTO sHayeHnä. Tax, HaNp., mPM TayccoBoM 
ABymapamerpnyeckoM pacıpezeleHun OH CONePp:KUT BOBpaCcTamımme CTemeHn 7” OÖEIYHOTO 
KOPpeAIHOHHOTO KoahhuuueHutar. C APyrofi CTOPOHBI, 1A 11000r0 3aNaHHOTO CHEKTPA Pac- 
mpemeneHunä YACTOT He HMeeTcH, & KOPPelAmmOHHBIe BHAYEHHA MOISKHEL ÖBIT OTPAHHYEHE 
Ha H3BECTHEIE ODNACTM. 


Um den Korrelationszusammenhang zwischen zwei Größen x und % festzustellen, für 
welche eine zweiparametrige Verteilung w(x,y) besteht, dienen Ausdrücke der Form 


ID Elena eier. 2. ll). 


der dann erhalten wird, wenn für f(x) und g(y) die beiden linearen Funktionen 
c—a | y—b 
le : 
verwendet werden, wobei die Konstanten a und 5 die beiden Mittelwerte und die Konstanten s 
und i die beiden Streuungen der Verteilungen 
e | u(2)= [w(a,y)dy ° bzw. a) eny)du 2... () 
bedeuten. Auch die zwei Pearsonschen Korrelationsverhältnisse fallen unter die genannten 
 - Größen, nämlich das eine für 
j | 2” —4 
Ze 


u; 


w(z,y) da, 


%— 4 
S 


ud  gW)= | 


Ns 


Die bekannteste und einfachste Größe dieser Art ist der gewöhnliche Korrelationskoeffizient r, 


und er N Ro | 


UR 
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das andere für r i BR: 
ar Ki VET ‚Yy)dy. 
W= und Soda) 7 vandy 
Die Frage nach einem Korrelationsmaß, welches den Zusammenhang zwischen den Größen & 
und y so deutlich als möglich anzeigt, führte den Verfasser dazu, das Variationsproblem 
K= [[f(e):g(y)-w(w,y)dedy=Maximum ..... eh, 
zu studieren!). Hierbei wurde von den unbekannten Funktionen (x) und g(y) verlangt, daß 
sie die Bedingungen 
[f(@)w(e) de — 0, Sswwly)dy=0 .:......- (958 
und i 
[S’(o)w(e) de=1, [eW)ws(y) Warner ta 


erfüllen, wie dies z. B. die linearen Funktionen (2) tun. Die Lösung dieses Variationsproblems 
ist die sog. „‚Maximalkorrelation‘, ein Korrelationsmaß, das u.a. den Vorzug besitzt, invariant 
gegenüber Merkmalstransformationen zu sein. ; 

Es gibt allerdings nicht nur ein einziges Maximum der Größen nach (4) sondern eine ganze 
Folge von Extremwerten. Sowohl für f(x) als auch für g(y) gilt:nämlich eine homogene Fred- 
holmsche Integralgleichung, deren Kern mit der Verteilungsfunktion w(z,y) zusammenhängt. 
Die Eigenwerte dieser beiden Integralgleichungen sind dieselben und mit ihrer Hilfe lassen sich 
deren Eigenfunktionen paarweise einander zuordnen. Ist A, der n-te Eigenwert und sind f, (=) 
bzw. g9,(y) die zugehörigen Eigenfunktionen der beiden Integralgleichungen, so nimmt mit ihnen 


1 ; a 
das Integral (4) einen Extremwert an, der K,— — beträgt. Wir vermerken noch für nach- 


folgende Anwendung die Formeln, welche die genannte Zuordnung der beiderlei Funktionen 
leisten. Sie lauten 


Inu y)dy= KW); [row y)de=Kng.(y)urly) -» - - (N). 


Da es im allgemeinen eine abzählbar unendliche Menge von Eigenwerten A, gibt, von 
denen auch mehrere miteinander übereinstimmen können, gibt es also auch eine Mannigfaltigkeit 
von extremen Korrelationswerten K,„. Sie bilden das ‚‚Korrelationsspektrum‘“‘ der vorliegenden 
Verteilung w(z,y). In der angegebenen Arbeit wurde vom Verfasser vorgeschlagen die engste 

dieser Korrelationen, die demnach zum kleinsten Eigenwert gehört, als „Maximalkorrelation‘ K 
> schlechthin zur Charakterisierung des in Frage stehenden Korrelationszusammenhangs zu 
verwenden. 

Die Brauchbarkeit dieses Vorschlags wurde später von Friede und Münzer bestritten ?). 
Diese Autoren betrachteten Beispiele für Korrelationstabellen, bei denen einerseits aus Sym- 
metriegründen der gewöhnliche Korrelationskoeffizient r sich zu Null ergibt, anderseits aber 
zahlreiche Felder in bestimmter Anordnung leer sind, wodurch die Maximalkorrelation X den 
Wert Eins annimmt. Wir werden im folgenden noch mehr derartige Beispiele antreffen, deren 
Besonderheit darin besteht. daß der kleinste Eigenwert, d.h. die engste Korrelation nicht zu den 
einfachsten, monotonen Eigenfunktionen gehört, sondern zu höheren Eigenfunktionen. Man 
kann jedoch nicht erwarten, daß es möglich ist, alle derartigen Besonderheiten mit Hilfe einer 
einzigen statistischen Maßzahl zu erfassen. Vielmehr liegt es nahe, einen Schritt weiter zu 
gehen und zu diesem Zwecke das ganze Korrelationsspektrum zu betrachten. 


Im folgenden werden wir einerseits für die wichtige zweiparametrige Gaußsche Verteilung 
das Korrelationsspektrum kennen lernen und an Hand dieses Beispiels Gründe erfahren, denen 
zufolge die Korrelationen höherer Ordnungszahl recht allgemein rasch an praktischer Wichtigkeit 
_ verlieren dürften, so daß die oben definierte Maximalkorrelation doch von überragender Be- 
deutung wäre. Andererseits wird uns die Aufgabe beschäftigen, zu vorgegebenen Korrelations- 
spektren die möglichen Verteilungen kennen zu lernen und zu erfahren, ob Korrelationsspektren 
beliebig vorgeschrieben werden können oder gewissen Einschränkungen unterliegen. Es ist 
jedoch nicht daran zu denken, das letztgenannte Problem hier in großer Allgemeinheit zu 
behandeln; wir werden uns vielmehr auf die Betrachtung einiger einfacher aber typischer 
Fälle beschränken. 
. Wir beschließen diese einleitenden Betrachtungen mit der Aufstellung einer für das Folgende 
wichtigen Beziehung. Gehen wir von der Annahme aus, wir besäßen für die noch unbekannte 
Verteilung w(z,y) die beiden Funktionen w,(x) und w,(y) nach Gl. (3). Außerdem seien das 


1) Z. angew. Math. Mech. 21 (1941), 8. 364—379, 


or ?) Z.angew. Math. Mech. 28 (1948), 8. 158—160, siehe auch H, Richter: Z. angew. Math. Mech, 29 (1949), 
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Korrelationsspektrum und die beiden Systeme von Eigenfunktionen f,„(x) und g„(y) bekannt. 
Die letztgenannten Funktionen erfüllen die Orthogonalitäts- und Normierungsbedingungen 


SH.) le) z w(&) dx SR [0 wenn mn 


Es gibt dann für w(z,y) eine Darstellung der Art (vgl. a.a. O., S. 374), 


CH LION AO EEK IOH TOR 
4,9= 

wofür Multiplikation beider Seiten mit f„(2)9,(y) und Integration über alle & und y den Koeffi- 

zienten 6,,,„ liefert gemäß der Formel - 


mn = IA) : 9n(Y) ; w(, Y) da dy. 


Da aber die Funktionen f(x) und g(y) hier nicht ein beliebiges System von Orthogonal- 
funktionen sind, sondern die Eigenfunktionen, für welche die Beziehungen (7) gelten, lassen sich 
die Koeffizienten c,„,„ ausrechnen. Es ist 

0 fü 

nn [Inte) [0.: Wut. ayda— [In Kıhlduda=g, " »Ln 

Daher wird die in Rede stehende Entwicklung 


vwYy)=wm@)u AR) ay)+t}:--:--- (9). 


Nach dieser Formel wäre die Verteilung w(z,y) aufzubauen. wenn das Korrelationsspektrum 
und die beiden Systeme von Eigenfunktionen vorgeschrieben sind. Es kann jedoch der Fall 
eintreten. daß die rechte Seite von Gl. (9) negativ wird, während die Funktion w(x, y) als Wahr- 
scheinlichkeitsverteilung nur positive Werte annehmen kann. Durch diese Bedingung wird bei - 
vorgegebenen Eigenfunktionen der Spielraum der Größen K, eingeengt, so daß das Korrelations- 
spektrum nicht ganz beliebig gewählt werden kann. Wir werden unten für diese Tatsache eine 
Reihe von Beispielen kennen lernen. 


Mittels der Beziehung (9) kann man auch sofort angeben, welchen Wert das Doppel- 
integral (1) annimmt, wenn andere als die Eigenfunktionen für f(x) und g(y) eingesetzt werden, 
Mit Rücksicht auf die Bedingungen (5) und (6) ist für diese beiden Funktionen folgender Ansatz 


zu machen: 
So)=ahA@W)tah@d+:: mt ++ =L| 


Ba ER 10). 
WENN): mit Hit — | an 

Hiermit wird das Doppelintegral (1) | 
y== SS) -ıy) way)dedy=abK+mbK,; +"... >. (11). 


Ist nun K=Max(K;) der größte der Korrelationswerte, d.h. die Maximalkorrelation, so folgt 
mit Rücksicht auf die Schwarzsche Ungleichung sofort für die Größe r nach (11) folgende 


. Ungleichung 


r=3ab,K, <Kiab, <Re) No yo—K,. 
i—1 i=1 1 1 


Werden also andere als die zur Maximalkorrelation gehörigen Eigenfunktionen als f(x) und g9(y) 
zur Berechnung eines Korrelationskoeffizienten gemäß Gl. (1) verwendet, so fällt das Ergebnis 
stets kleiner als K aus. 

Es sei noch vermerkt, daß alle in diesem Abschnitt angegebenen Beziehungen in leichter 
Abwandlung auch für Korrelationstabellen mit % Spalten und I Zeilen gelten. Diese endliche 
Theorie ist im Werk des Verfassers „‚Zahl und Wirklichkeit“, 2. Aufl., Heidelberg 1949, S. 366ff. 


£ 


ausführlich dargestellt. ; 


Das Korrelationsspektrum der zweiparametrigen Normalverteilung 
Unter den Wahrscheinlichkeitsverteilungen für zwei Veränderliche spielt die wichtigste 


BE Rolle die zweiparametrige Gaußsche Normalverteilung.-Es ist daher von Interesse zu erfahren, 
was die hier vorgetragene Theorie für dieses Beispiel zu sagen hat. 


Wir können von der normierten Form dieser Normalverteilung Gebrauch machen, bei der 


E = 2 2 die Erwartungswerte a und b gleich Null und die Streuungen s und i gleich Eins sind. Der einzige 


12 


noch frei verfügbare Parameter ist dann bekanntlich der Korrelationskoeffizient r und die 
Verteilung lautet 


Sf. ®—2r2y+Y ra 
h ’ = ER ex st BETT TE EHEE NN LT. . . * . . 12 
(z ) 21 —r 2 21 —r) =) | = ). = 
Be: _ Integration dieser Funktion über y bzw. & nach der Vorschrift (3) liefert = . E r= = x; =. A 
nn 1 \ und Ay)= DE or er 
® % —_ ex ——) ın eg a NE 
TE \ be y2ar | s y y2 T _ | = h es En 


Es soll nun gezeigt werden, daß die Orthogonalpolynome steigenden Grades für ‚diese i 
Belegungen, das sind im wesentlichen die Hermiteschen Polynome, die Eigenfunktionen des” - 
Korrelationsproblems für die Verteilung ur darstellen- Für diese normierten Polynome eilt, Y 


x 5 ® die Formel WE 
Ki 1" ae\ .d" Bee N 
9 (2) = Ti exp 5 BR ein a Ene „9. we: 
Ohne den Normierungsfaktor Er geschrieben lauten sie Pan RR re 
n! ee 


2a) ar ee + He Vene n. (ie), 


wobei diese Reihe bei k 7: ; oder Et ur abbricht, je nachdem n gerade oder ungerade ist. 

Die ersten dieser Polynome sind - 
"p=l, H=, HB=eR—l, B=R—35 mer—6R+3, HB=R— 10a + 1de. 

Für nachfolgenden Gebrauch vermerken wir noch die Formel für die m-te Ableitung von Pr (2): | 


px) ee [n]m gn-m _ DE ER + — ..t (— 1)* ee Kl an-2k—m E- (14b). 


Um zu bestätigen, daß die Polynome p, die Eigenfunktionen ler Verteilung h(z, y) sind, 
muß das Kine et 


Bi: = * | Lm@) hy) dz 
berechnet werden. Zu diesem Zwecke wird zuerst eine Formel für 
SE (2 — ry)* -h(z,y) dx 


Be bereitgestellt, Man kann den Exponenten von h(z,) folgendermaßen umformen: 


R—2rayty? _ (® rer) (._ ry)? 
| 2(1—r?) a 2(1-—r?) 2(1—r) 
Damit wird ! a 
2 de 
Beyynanger lem resp (— Gh) as 
re RE Da 
Pr ) Pal- Fe Er) er 


Nun ist, nach Taylor um die Stelle 2&—=ry iirka, 


en le ee 


- 
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und daher 


oo 


?n(&) h(z, y) de = 


= [m 


Es handelt sich jetzt nur noch darum, den Klammerinhalt rechts auszurechnen, der, wie man 
sofort sieht, ein Polynom n-ten Grades in yist. Der Koeffizient von y" lautet wegen (14a) r”; 
derjenige von y*-? ist mit Rücksicht auf (14b) 


ar it 


tee met, 


derjenige von En Feist 


In 


— (1—r2) —— _— [n]* mtl — Peg E 


22.21 DT ARSCLNE 
= ee a A142, Pb; 
\ ‘- und allgemein ist der Koeffizient von A 
| Amer 
s N a ne 


= (pP Bet rl) (1-1) +) (r2)2 .. ) ae er ir 


"Das Ergebnis wird hiermit 


| la da alu) (y— a u 1100.20} 


Vergleich mit Gl. (7) zeigt nun, daß in der Tat die Funktionen p,(x) und p,(y) zusammen- 
gehörige Eigenfunktionen des in Rede stehenden Korrelationsproblems sind und daß der zu- 
gehörige Korrelationskoeffizient n-ter Ordnung K,„=r” ist. Das Korrelationsspektrum der 
Gaußschen Verteilung ist demnach besonders einfach, denn die aufeinander folgenden Korre- 
lationskoeffizienten sind nichts anderes als die steigenden Potenzen des gewöhnlichen Korre- 
lationskoeffizienten r, der hier zugleich die Maximalkorrelation angibt. Daher ist man hier auch 
wohl berechtigt, den ersten Korrelationskoeffizienten r als charakteristisch zu betrachten. 


> E. Erinnert man sich an die normierten Polynome (14), so gewinnt man sofort für die Normal- 
verteilung h(z,y) noch die bemerkenswerte Beziehung 


3 2 0 für m £n 

E Im: mW Hayden mumın ee a5), 
3 _ sowie auf Grund von Gl. (9) die Reihenentwicklung 

E: X h(z,Y) _— 

Er S Ka) I Ta@a)+z na) +: 3 ROYRO) ze - (16). 
E 5 Soweit der theoretische Sachverhalt über das Korrelationsspektrum einer normalen Ver- 
zn teilung. Um nun auch ein Urteil über die praktische Bedeutung der höheren Korrelations- 


koeffizienten zu gewinnen, wenden wir uns nun der Frage zu, wie sehr diese die Neigung haben, 
sich bei einer Stichprobe zu reproduzieren. Die Antwort ist in folgendem Satz über die Streuung 
der Maximalkorrelation enthalten:°) Wird in einer Stichprobe mit m Elementen eine Verteilung 
= angetroffen, für welche die Maximalkorrelation K beträgt, so ist in der sehr viel größeren Gesamt- 
- heit eine Maximalkorrelation zu vermuten, deren Erwartungswert näherungsweise mit K über- 
ee einstimmt und deren es en en wird durch 
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wobei f=/(x) und g=g(y) diejenigen Eigenfunktionen sind, welche zur Berechnung des extremen 
Korrelationskoeffizienten K führen. £ 

Zum vorliegenden Beispiel ist Gl. (17) für die Verteilung h(®; y) nach Gl. (12) und für 
K,=r", [=p,(x), g=p„(y) auszuwerten. Dies geschieht am einfachsten, indem man das 
Quadrat der Klammer unter dem Integral in einem Ausdruck der Form PX} Pl) ver- 
wandelt und darauf unter Berücksichtigung von Gl. (15) gliedweise integriert, wozu übrigens 
nur die zahlenmäßige Kenntnis der Koeffizienten c;; erforderlich ist. Das Ergebnis der etwas 
langwierigen Rechnung für n=1, 2 und 3 lautet: 


für Kh=r: Ym-Str(K)=1-—1r? 
für R=r: Ym-Str(K)=2(1—r)yl+10r +37 STIER 
für R,—=r: Ym-Str(K,)=6(1—r))1 +20r +93 +47 r°+ 107° 


In Bild 1 ist der Verlauf der mit Ym multiplizierten Streuungen dieser Maximalkorrelationen 
erster bis dritter Ordnung über K aufgetra- 
20 gen. Das auffallendste Ergebnis ist die rasche 
Ym:Str (K) Zunahme der Unsicherheit bei den Korrela- 
tionen steigender Ordnung. Es haben daher 
etwaige Abweichungen der höheren Korrela- 
tionen von den theoretischen Werten r” nicht 
viel zu bedeuten. Das weitaus größte Gewicht 
hat doch normalerweise die Maximalkorre- 
lation K=K,, das ist hier bei der Gaußschen 
Verteilung der gewöhnliche Korrelationskoeffi- 
zient r. 

Da überdies auch die logarithmischen 
Normalverteilungen für zwei Parameter das- 
selbe Korrelationsspektrum besitzen, ist wohl 
anzunehmen, daß bei den weitaus meisten in der 
Praxis vorkommenden Verteilungen die Korre- 
lationen steigender Ordnung immer lockerer 
werden, indem ihre Koeffizienten ungefähr 
nach dem Potenzgesetz r* abnehmen. Man 
wird dieses Verhalten als ‚‚normal‘ betrachten 

rer dürfen, während der gegenteilige Fall, wo Korre- 
Bild1. Streuungen der Größen K,, K, und K, lationen höherer Ordnung besonders groß ausfal- 
len, als ausgesprochene Abnormitätgeltenkann. 


Korrelationstabellen zu vorgeschriebenem Korrelationsspektrum 
Wir wollen nun umgekehrt vorgehen und Wahrscheinlichkeitsverteilungen herstellen, die 
neben anderen Bestimmungsstücken ein vorgeschriebenes Korrelationsspektrum besitzen. Dabei 
beschäftigen wir uns im folgenden nicht mehr mit stetigen zweiparametrigen Verteilungen w(z, y), 
sondern beschränken uns auf die Betrachtung einfacher Korrelationstabellen mit wenigen 
Spalten und Zeilen. 
Einleitend sei vermerkt, daß es bei einer Korrelationstabelle mit % Spalten (x-Werte) 
‚und I Zeilen (y-Werte) k—1 nicht konstante Orthogonalfunktionen f}(&) bis fz-ı (x) und I—1 
nicht konstante Orthogonalfunktionen g, (2) bis 9-ı(x) gibt. Ist k<I, so‘können höchstens 
k—1 Korrelationswerte K, von Null verschieden sein. (Vgl. „Zahl und Wirklichkeit“, S. 376.) 
1. Bernoullische Korrelation: Der einfachste Fall ist der einer Korrelationstabelle mit 
‚nur zwei Spalten und zwei Zeilen; er wird bekanntlich als Bernoullische Korrelation bezeichnet 
Ganz allgemein lautet die Korrelationstabelle 27 


An Raı | Qı 
a Rıs LEN an 
Pı Pa 1 


Die Matrix der vier Häufigkeiten Ah;;, die zusammen die Summe 1 ergeben, ist z.B. bestimmt, 
wenn die Größen 7,=?, q,—=qg und Ah, vorgegeben sind. Es gibt also 008 derartige Matrizen. 
Das System der f(x) reduziert sich hier auf eine einzige Funktion, von der auch nur zwei Werte fı 
und f, interessieren. Für f, und f, gelten entsprechend den Beziehungen (5) und (6) die beiden 


Gleichungen: 
| | Amt+hm=0; Sat+tfm=t; 


j 


\ 
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deren Lösungen, mit 9 —=p und 9—=1—p geschrieben, lauten 


und entsprechend ist 


H= IE y=— ee ee ee (19a). 


Die Reihenentwicklung nach (9) ist hier, da es nur einen von Null verschiedenen Korrelations- 
koeffizienten r gibt, einfach 
;=pgyltrfig), 


oder ausgeschrieben mit der Abkürzung 


N EN En A (20) 
| nasse ee 
ee (21). 


In dieser Darstellung enthält die Matrix die drei Parameter p, gund r. Während p und g beliebige 


“positive echte Brüche sein dürfen, kann r nicht beliebig zwischen —1 und +1 vorgegeben werden, 


denn die vier Größen der Matrix (21) müssen als Wahrscheinlichkeiten alle positiv sein. Dies 
tritt nur ein, solange r größer als die größere der beiden Zahlen 


-/ pq nt peg 
1—p)1—g En Y vq 


und kleiner als die kleinere der beiden Zahlen 
pl—q d—p)gq 
ne und Ve 
1079 % p(l—q) 


Die obere Grenze für r liegt bei +1, wenn p9=gq ist, und die untere Grenze liegt bei —1, 
wenn p+g=1 ist. Nur im Falle 


ist. 


1 
kann r alle Werte zwischen —1 und —1 annehmen. Die Matrix reduziert sich hierfür auf 
a 
4\11—r. 1-#r 
Die obere und die untere Grenze für r liegen $ymmetrisch, wenn p =, oder q =; ist 
1 1 
Sie betragen z. B. für = 7. und gs 7 
el VE N 
1—q a 1—q 
Sind p und q von 5 verschieden, so ist der Bereich für r unsymmetrisch. So lauten z.B. für 
»=0,3 und q= 0,4 die Bedingungen | | 
0,3 -0,4 0:3.°0,6%7 
r> Vor 058: Me 0,7 GG 


2. Beispiele für Matrizen mit drei Reihen und Spalten: Bei der Betrachtung drei- 
reihiger Matrizen beschränken wir uns auf den Fall 


1 1 
m—cont=-7; g—wnt=-7, 


ee da uns mit den bekannten Orthogonalpolynomen für diese Belegung sofort ein geeigneter Satz 


von Eigenfunktionen f,(x) bzw. 9,(y) zur Verfügung steht. Dies bedeutet eine wesentliche 
Einschränkung der Allgemeinheit, denn während die allgemeine Korrelationstabelle mit je drei 


Be Reihen und Spalten acht willkürliche Parameter enthält, wird durch die Vorgabe der p; und q; 
- über vier Parameter verfügt, so daß die noch möglichen Lösungen der gestellten Aufgabe eine 
- vierparametrige Mannigfaltigkeit darstellen. : 
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je i nichtkonstanten Funktionen, 
Das System der f„(x) bzw. g,(y) besteht hier nur aus zwei nie 
von denen überdies nur drei Werte interessieren. Diese sind in folgender Übersicht zusammen- 


gestellt. 3 3: 
ha=-Vz x +Ya a RENEN Be 


na=+)/g 12 +5 uw 


Für diejenigen Korrelationstabellen, für welche dies die Eigenfunktionen des Korrelations- 
problems sind, gilt nach (9) die Entwicklung 
1 


h(&;, y) = 9 (1 + Kıfı(&)g(y) + Beh(®) 9(4;)) 5 
oder ausgeschrieben 
z 3 1 3 1 K 
[145 Kı+z&: 1—KR, 1-ZKırz \ 3 | 
!| 1-K, 1+2K, 1—R, ER (23)... 0 
. 3 1 3 1 % . 
1—-zKırz&: 1—K, 1+5KırzR, 


In dieses Ergebnis sind noch die beiden willkürlichen Parameter K, und K, eingegangen. Es 
Fa stellt auch für unser spezielleres Problem noch nicht die allgemeine Lösung dar; wir wollen es _ 
B: aber zunächst diskutieren. | 
NER „2 Da alle Elemente der Matrix (23) Häufigkeiten darstel- 


len und daher keine negativen Werte annehmen können, 
unterliegen die beiden Parameter X, und K, gewissen Ein- 
schränkungen noch über die Bedingung hinaus, daß sie als 
Korrelationskoeffizienten die Werte +1 nicht überschreiten 
dürfen. Der Bereich für X, und K,, welcher durchweg 
K, positive Häufigkeiten in (23) zur Folge hat, ist das schraf- 
fierte Trapez von Bild 2. In diese Zeichnung findet sich auch 
die Parabel K,—=K? eingetragen. Diese Kurve liegt im 
Inneren des erlaubten Bereichs. Nach dem obigen Ergebnis 
BER, nnz für die Gaußsche Verteilung dürfte es als ‚„‚normal‘ gelten 
können, wenn K, und K, wenigstens angenähert in dieser 
Bild 2. Bereich für K, und K, für die Matrix Beziehung zueinander stehen. Auffallend ist, daß für R,—=0- 
Ka a die Größe K, nicht mehr beliebig gewählt werden kann, son- 


dern auf den Bereich zwischen— und + = beschränkt bleibt. Dagegen kann K, alle Werte 


- zwischen —1 und +1 annehmen, wenn R,—1 ist. - In diesem interessanten Fall reduziert sich - " 
- die Matrix (23) auf 


A OL GEN.c, | ER > Be 
it DR en a DI RR a 
Dies ist ein besonders einfaches Beispiel der von Friede und Münzner (s. oben) kritisierten Art. 
Obwohl hier die Maximalkorrelation den Wert Eins aufweist, wird für K,—0 die gewöhnliche: ‘= 
Korrelation Null. In Wirklichkeit sind jedoch die beiden Veränderlichen keineswegs voneinander Er 
_ unabhängig. Vielmehr liegt in der Tatsache der vier leeren Felder eine sehr starke Bindung, de - 
zur Folge hat, daß durch zwei Faltungen sogar die Korrelation K,—1 zutage tritt. Um solche 


abnorme Korrelationstabellen zu analysieren, ist es allerdings erforderlich, das vollständige 
_ Korrelationsspektrum zu betrachten. Br . Pe ng 


. Zum Abschluß dieses Beispiels von Gl. (23) betrachten wir noch den Fall KR SE 
Damit’reduziert sich die Matrix auf SE RR DR EN. 2 | 
i+2r I—r 1-—r HR 5 Be re Be 
: 3.1.1 er er ae rd Fi « (23b). Fa 
rs: Ir 18 bar) 
“ Dies ist, die naheliegende Verallgemeinerung des Ergebnisses (21 a). Hier tritt für r=0 voll- Se 


er, ? De " r 
N - 
i Ser 


NN 
\ y . 
. 


EN 
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ständige Unabhängigkeit und für r—=1 vollständige Abhängigkeit ein. Der niedrigste Wert, der 
für r möglich ist, beträgt = 

Wie schon oben festgestellt worden ist, liegt mit der Matrix (23) noch nicht die allgemeinste 
Lösung der gestellten Aufgabe vor. Die Einschränkung besteht in der Annahme, daß gerade die 
Funktionen (22) zu Eigenfunktionen werden sollen. Es gibt aber neben (22) noch andere Paare 
von normierten Orthogonalfunktionen zur gleichen Belegung, die aus ihnen durch geeignete Linear- 
kombinationen zu gewinnen sind. Diese lauten 


F,=fic0sa+fssin o; H=gcosß+ sin; | 


EL RR: (24). 
F,—=fısina —facos a; G,—gsinß—gcosß J 
Mit diesen Größen ist dann 
1 i 
h(&,,y)= 9 (1 + KıF,(&) 61 (y)+ Kr Fr(e;) G,(Y;)) SI (25) 


die gesuchte allgemeine Verteilung mit den vier erforderlichen Parametern K, K, & und Pf. 


Zwei Beispiele mögen einen Begriff von der Mannigfaltigkeit der in diesem Ergebnis 
liegenden Möglichkeiten geben. 


a b—45: 


Die Korrelationstabelle lautet: 
1+ 0,134 K,+ 1,866 K, 10,366 K,— 1,366K, 1—0,5K,—0,5K, 
5 1-+ 0,366 K,— 1,366 X, 1+ K, +, 1 — 1,366 K, + 0,366 K, 
1—0,5K,—0,5K, . 1-—1,366 K, + 0,366 K, 1-+ 1,866 X, + 0,134 K, 


Bild 3 zeigt für diesen Fall den Bereich der möglichen Parameterwerte X, und K,. Auch hier 
ist dieser Bereich ein Viereck, das jedoch gegenüber Bild 2 verschoben ist. 


«= 0°, B= 90°: 


Während bei den bisherigen Beispielen die Funktionen f„(x) und g,(y) dieselbe Reihenfolge auf- 
weisen, ist hier #},—=@, und F,—=@;, wodurch eine merkwürdige Störung der Symmetrie hervor- 
gerufen wird. Die in Rede stehende Matrix ist 


rent 1+J3 K, een, 

1 # 

g|jitY3R, 1 1—Y3 &, 
ee ee 


K, 


+1 


i Bild 4. Bereich für K, und K, 
a m im Falle & = 0°, ß = 90° 


N 


Den Bereich der möglichen Werte K, und K, zeigt Bild 4. Hier begegnet uns zum erstenmal 
eine Figur, die an keiner Stelle den Rand des gestrichelten Quadrats erreicht. 


"Es verdient darauf hingewiesen zu werden, daß bei dem letzten Beispiel («—0), Bß=%°) 


; ‚stets der gewöhnliche Korrelationskoeffizient r—=0 ist. Jedoch kann von stochastischer Unab- 


2 


N 
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hängigkeit keine Rede sein, wenn Korrelationstabellen so merkwürdige Besetzung zeigen wie 
die folgenden: 


1 rl 
— KR, —= — = 0,577: Kıe—ı  ) 
K, K; 3 1 y3 s 2 
70 .2..1 1. 223 
- 22,0 as 4.2 -0.5 
9 £ 18 
ID. 2 12 
3, Matrizen mit vier Reihen und Spalten: Auch hier beschränken wir uns aus dem- 
selben Grunde wie bei den dreireihigen Matrizen auf den Fall 
1 1 
| p;— const —a g; = const = 7: 
Im übrigen können wir uns nach der ausführlichen Besprechung des vorangehenden. Beispiels 
b- kurz fassen, indem wir die Orthogonalpolynome steigenden Grades zu diesen Belegungen unver- 
= ändert in ihrer natürlichen Reihenfolge als Eigenfunktionen verwenden. Es ist dann 
er =3 1720 432 
= j 1 ee — a — ; 
= £ Fı( ) 5 y5 y5 y5 9(Y;) 
& a)=—ı + Tea (26). 
ae .. 


Die Reihenentwicklung nach (9) lautet hier 
h(2,9;) = = (1 + K, f(&)9ı (4) + Ka fe(&) 9(Y) + Ks Fs(&) 9s)) 


oder ausgeschrieben 


(eisen Ks; 543 KR BR, SR, 5SK SR SR 5 OR ıSR 


ve = 9(9;) 


1 


54H5K,SK,3K,, 54 KA5K, IK, 5 Kıi5R, IK, SaK ES 
80 


53 KB, 5 Ry43K;,"5— Kıl5R, OK 5 RK SK OR Ba an 
5-9K,+5Ry- K,, 5-83K,-5K,t+3K,, 543K,5K,3K,, 549K,4+5K,t K, 
Für den Fall K,=0 wird der Bereich der. möglichen Wertepaare K,, K, durch Bild 5 
angegeben. Interessante Grenzfälle sind die obere und die untere Spitze der drachenförmigen 
schraffierten Figur. Die zugehörigen Matrizen sind von dem 
abnormen Typus, der bereits oben besprochen worden ist. 


K=-K=0; B=+1; B=K,=0; B=—1; 


A a BE 
110 Ber 1.1 120% 00 
| TE TEE Be 

1: Soest 04.1 Ve 


Allgemein für X, 0 ist der Verträglichkeitsbereich für 
K,, K, und K, die in Bild 6 dargestellte räumliche Figur. Als 
Bild 5. Bereich für X, = 0 Schnitt mit der (K,, K,)-Ebene ist die drachenförmige Figur von 
Bild 5 wieder zu erkennen. Auch die räumliche Kurve KK: 
K;—K}, die den Verhältnissen des Normalfalls entspricht, ist eingezeichnet. Sie verläuft auch 
hier im Innern des dargestellten Körpers. 


Interessante, abnorme Korrelationstabellen gehören insbesondere zu den Punkten der 
oberen Kante K,=1, KR=K,=r. Dafür ist die Matrix 


i+r- Teer 
0 N > LE 
810 sen 
1—r 0 0 .1-+r 


ee : ae an ihr ist, daß beide Diagonalen mit verschiedenen, konstanten Werten besetzt 
Br "sind. Pe re 


Bd El 


> 


TB 


ET, Eingegangen, den 14. März 1951. 
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Ein weiterer, interessänter Fall liegt endlich vor, wenn mit Kı= K,=K,=r alle vier 
Extremalkorrelationen übereinstimmen. Das Ergebnis 


1-+3r 1—r 1—r 1—r \ 
a; 1+3r 1—r 1—r 
16| 1—r 1—r 1+3r 1—r | 

1—r 1—r 1—r ee 


ist die einfache, zu erwartende Verallgemeinerung von (21a) und (23b). 

4.Bemerkung über allgemeine 
Korrelationstabellen: Es wurde oben 
am Beispiel der Korrelationstabelle mit je 
drei Spalten und Zeilen gezeigt, wie die 
acht Parameter zustandekommen, die das 
Problem in seiner vollen Allgemeinheit be- 
stimmen. Diese Überlegung sollnun noch 
auf den Falleiner Korrelationstabelle mit 
k Spalten und } Zeilen verallgemeinert 
werden, wobei wir ohne Einschränkung 
der Allgemeinheit annehmen können, daß 
k<I sei. 

Zunächst sind (k— 1)vonden %k Wer- 
ten p; und (—1) von den l Werten q; ver- 
fügbar. Weiter gibt es (k—1) möglicher- 
weise von Null verschiedene Korrelations- 
koeffizienten K,. So weit sind also 
(2k-++-1—3) Parameter festgestellt. Nun 
gibt es zur Belegung p; ein System von 
(E—1) nicht konstanten, normierten 246. Mäunlicher Basic I Ai, K, uad Fr be der Korrelation 
Orthogonalpolynomen für x, zur Bele- 
gung g, aber ein solches von (I—1)nicht konstanten, normierten Orthogonalpolynomen für y. Aus 
diesen sind die beiden Serien von je (k—1) Eigenfunktionen f,(x) und g,(y) aufzubauen. 

Da die Eigenfunktionen f,(x) zueinander orthogonal sind und normiert sein sollen, und 
da für sie (k—1) Basisfunktionen zur Verfügung stehen, können (k—2) Konstante bei der Auf- 
stellung von f(x) gewählt werden, (k—3) Konstante bei der Aufstellung von f(x) usw. Die 
Gesamtzahl der für die Bildung der (k—1). Funktionen verfügbaren Konstanten ist daher 


+94: +@-H=&-1(3-1). 


Entsprechend können (I—2) Konstante bei der Aufstellung: von g,(y) gewählt werden, 
da es hierfür (—1) Basisfunktionen gibt, (—3) Konstante bei der Aufstellung von 9,(y) usw. 
Die Gesamtzahl der verfügbaren Konstanten für die Bildung der (k—1) Eigenfunktionen ist 


ÜyHE- Ir. +0 y=0-Yli-3—ı). 


Nun ist schließlich die Gesamtzahl der genannten Parameter 
k k 
@2+1-9+6@-1(3-1) +@=yl1-3-1) Sa 


Das ist genau die Zahl der Häufigkeiten h,; vermindert um 1 infolge der Bedingung, daß die 
Summe dieser Häufigkeiten Eins ergeben muß. De 
Schlußbemerkung: Das Korrelationsverhalten einer zweiparametrigen Häufigkeits- 
verteilung läßt sich grundsätzlich nicht durch Angabe eines einzigen, irgendwie gearteten Korre- 
lationskoeffizienten beschreiben. Vielmehr ist zu diesem Zwecke prinzipiell das Studium des 
zur Verteilung gehörigen Korrelationsspektrums erforderlich, das bei einer Korrelations- 
tabelle mit % Zeilen und Spalten aus (k—1) durch Extremumseigenschaften ausgezeichneten 
Korrelationswerten K, besteht. Trotzdem dürfte die Kenntnis des ersten bzw. größten dieser 
Werte für die Praxis im allgemeinen ausreichen, denn erstens sind bei den besonders häufig vor- 
kommenden Normalverteilungen die Größen K, nichts anderes als die aufeinanderfolgenden 
Potenzen des ersten Wertes K, und zweitens werden die Korrelationen höherer Ordnung bald 
sehr unsicher, da die Streuungen der K;, mit wachsendem n stark zunehmen. Das Studium von 


(11,1) 


- Korrelationstabellen im Zusammenhang mit dem Korrelationsspektrum eröffnet einen Weg, um 


die mannigfaltigen Möglichkeiten zu überblicken und Abnormitäten besser zu verstehen. 
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KLEINE MITTEILUNGEN 


Die Zurückführung von speziellen linearen 
Integrodifferentialgleichungen auf gewöhn- 
liche Integralgleichungen. 

A. Es handelt sich um die Auflösung der folgenden 
linearen Integrodifferentialgleichung: 


b 
ya) -ASCRHNWHldE=fe). - WM. 


Dabei sollen bedeuten: 
1. N [y(x)] einen linearen homogenen Differential- 
ausdruck n-ter Ordnung in y(x) von der Form: 


N v 
N ua] = 2 te) (1a) 


d 
mit in <a, b) stetigen (reellen) Koeffizientenfunktionen 
(2), vd =0,1,2, *,n. 
2.@(x,£) einina <x, € Sbstetiger (reeller) Kern, 
dessen partielle Ableitungen nach x bis zur n-ten 
Ordnung einschließlich 


0a,E) PAS), 


Or 02° I 


ina<x,&<<b existieren und dort stetig sein sollen. 


3. f(x) eine in <a,b)> n-mal stetig differenzierbare, 
sonst beliebige (reelle) Funktion. 

B. Ein Lösungsweg besteht darin, die Gl.(1) durch 
n-malige Differentiation in ein Simultansystem von 
Fredholmschen Integralgleichungen für die Funk- 
tionen y(z), Y (2), y”(&),***, y®(x) überzuführen und 
dann das entstandene System in bekannter Weise 
auf eine einzige Integralgleichung zu reduzieren. 

Zweckmäßiger erweist sich aber ein Verfahren, wel- 
ches die Tatsache, daß in der Gl.(1) der Kernopera- 
tor @ und der Differentialoperator N gewissermaßen 
getrennt auftreten, ausnützt. Dieses Verfahren be- 
steht im folgenden: 


C. Infolge der Voraussetzungen 2. und 3. kann der 
Differentiator N auf die Integrodifferentialgleichung 
(1) angewandt werden, so daß sich ergibt: 


b 
Niy(@)]—A j De - 
a 


NG (x,E) 


9x 


—_ (2) 
"N [y(old& = NL[f(e)] 
Man setzt noch zur Abkürzung: 
NIyayl= Da) Ed mn; | 
yv=0 T 


Nlf@)]=Fe). 


Mit diesen Bezeichnungen Gl. (3) liegt in Gl.(2) die 
folgende lineare Fredholmsche Integralgleichung 2.Art 
für ®(x) vor: 


b 
D(e) — af H(w6) D()dE=Fla) .. (4). 


Die Integralgleichung (4) sei als die ‚‚begleitende 
Integralgleichung“ der Integrodifferentialgleichung (1) 
bezeichnet. Sie besitzt einen wegen der Voraussetzun- 
gen 1 und 2 stetigen Kern H(z,£), und auch F(x) ist 
nach Voraussetzung 3 in.<a, b) stetig. Der Ausartungs- 

fall H(x,8) = 0 sei ausgeschlossen. 

‚Die begleitende Integralgleichung a estattet nun 
die Auflösung der Integrodifferentialg eichung (1), 
d. h. mit Lösungen der Integralgleichung (4) sind die 
Lösungen der Integrodifferentialgleichung (1) durch 
eine explizite Formel gegeben. Zum Beweis unter- 


en wir den homogenen und den inhomogenen 
all, 


D. Homogener Fall: f(x) = 0, also nach Gl. (3) 
auch F(x) = 0. Die begleitende Integralgleichung ist 
dann ebenfalls homogen: 


b 
D(a)=Af Hia,E) D(e)dE. . . - (A). 


Es seien 
Ay Ass Ass 2 Aus wi (5) 


die Eigenwerte des Kerns H (x, £) der Integralgleichung 
(4), je nach ihrer Vielfachheit entsprechend oft auf- 
geführt. Die zugehörigen, der Eindeutigkeit halber 
über <a,b> normierten Eigenfunktionen seien mit 


D, (x), D;(x), D; (x), ee) Du), ii (6) 


bezeichnet. (Für einen r-fachen Eigenwert Au = Ausı 
= ++ = Ayyr-ı ordnet man jedem Au+, v% = 0,1, 2, 
+, r—1, je eine der r-normierten, voneinander linear 
unabhängigen Eigenfunktionen zu.) 

Es gilt dann der Satz |: 

Die Gesamtheit der Eigenwerte der homogenen 
Integrodifferentialgleichung (1) und die Gesamtheit 
der Eigenwerte Gl. (5) der begleitenden Integral- 
gleichung (4°) stimmen überein. . Ist ®,(x) die zum 
Eigenwert A, von H(z,£) gehörige Eigenfunktion, 
so ist 


ae Burn 


b 
Yule) Anl RE, Dule)dE ... M 


die zum Eigenwert A, gehörige, bis auf einen kon- 

stanten Faktor-bestimmte Eigenfunktion der homo- 

genen Integrodifferentialgleichung (1). 

Beweis: Essei y*(x) eine zum Eigenwert A* ge- 
Aotige Eigenfunktion der homogenen Gl.(1), es gelte 
also 


d 
he; a NLyrälas . . (1). 


Die Anwendung des Differentiators N auf diese Re- 
er ergibt unter Verwendung der Bezeichnungen 
1.(3): 


b 
N [y*(«)] RER H(&,$) N[y*(S)]lae. 


Da N[y*(x)] #0 (sonst wäre nach Gl.(1’) auch 
gegen die Voraussetzung y*(«)=0), muß N[y*(z)] 
die bis auf einen konstanten Faktor bestimmte, zum 
Eigenwert A* gehörige Eigenfunktion von H (x,£) sein. 
Die Eigenwerte der Integrodifferentialgleichung (1) 
sind somit auch solche des Kerns der begleitenden 
Integralgleichung (4°). 

Umgekehrt bedeute jetzt ®,„(x) die zum Eigenwert 
Au gehörige, normierte Eigenfunktion von H(x,£). 
Man bildet nach Gl.(7) die Funktion 


b 
Yulz) =Auf ClasE) Bule)dE nl) 


und zeigt, daß diese Funktion y„(x) eine Lös der 
Integrodifferentialgleichung (1) im homogenen Falle 
für A=), ist. Dies ist dann und nur dann der Fall, 
wenn die nach Gl.(7) erklärte Funktion y„(x) die 
Differentialgleichung ; 

Nly@)l= Due) ....., 8 
mit der stetigen rechten Seite ®,„(x) erfüllt. 


Nun ist aber unter Verwend Gl. (3), 
endm ung von (3) 


b 
N [yu(®)] MS H (2,8) Du($)dE = Dula) (8). 
Aus Gl.(7) und (8°) folgt: 
b 
Yu(z) — Auf @(z,E) N [yulOldE. 


(a) 
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Dabei ist y„(x)=0, da sonst nach GI. (8°) auch 
Du(z)=0-wäre, wohingegen ®,(x) doch eine Eigen- 
funktion von H(z,£) sein sollte. Da alle Eigenfunk- 
tionen der homogenen Integrodifferentialgleichung (1) 
die Form 61. (7) haben müssen, ist der Satz 1 in allen 
‚Stücken bewiesen. 

E. Inhomogener Fall: f(x) #0. 

Unterfall I: N[f(e)]=F(2) #0. 


Es sei A=/, kein Eigenwert von H(z,£). Es exi« 
stiert dann nach dem 1. Fredholmschen Theorem eine 
stetige, eindeutig bestimmte Lösung Ö,(x) der beglei- 
tenden Integralgleichung (4). Mit dieser Lösungs- 

funktion ®,(x) ergibt sich der Satz 2: 
Für A9=#Au (A. Eigenwert des Kerns H(z, &)) ist 
_ die Lösung der Integrodifferentialgleichung (1) ein- 

-  deutig gegeben durch 


b 
Yo(z) = As Se@5) Dle)ds+fle). . (9. 


- 


9: - Beweis: Man stellt wie beim homogenen Fall 
B ee fest, daß die durch G1.(9) definierte Funktion y,(x) 
Fe. der Differentialgleichung 


N[y(@)] = ©s(2) 

genügt. ZEREN € 
_ _ Die Einzigkeit der gefundenen Lösung yo (x) folgt 
' daraus, daß für alle Lösungen der Gl.(1) der zuge- 
hörige Differentialausdruck N[y(z)] die begleitende 
Integralgleichung (4) erfüllen muß. Da aber für 
AA. die begleitende Integralgleichung (4) nur eine 
_ eindeutig bestimmte Lösung besitzt, und die Lösungen 
von Gl.(1) notwendig die Form G1.(9) haben müssen, 
kann auch nur eine einzige Lösung von Gl.(l) exi-. 
 stieren. i Ay 
_ IstA=#,„ ein Eigenwert von H(x,£), so findet im 
Inhomogenen Falle F(x) #0 das 3. Fredholmsche 
- — Theorem auf die 61.(4) Anwendung. Den Rückschluß 
auf gegebenenfalls existierende Lösungen von Gl.(1). 
führt man analog wieoben. « 
Unterfall II: f(«) #0, F(o)=0. 
In diesem Falle setzt man 
ua) 2. =. (10) 


hal aus der inhomogenen GI. (1) die homogene 


ichung. 


Er .e a b er £ ee S 
u) SCmENTule]dE . 
z a Re 


ur 


a 


begleitende Integralgleichung von Gl.(11) hat 


 Dla)=A SH) BiHde... (M. 
Wr ee a a - en er Be 

rerte A „von H(x,£) wird nach Satz 1 
Integrodifferentialgleichung (11) ge- 


— 
. 


KOEHE. 


Pe: ee” : den Potenzen von n berechnen läßt. Zeichnet mansich 
danach die Kurve f(n) etwa für n=3+++12,5 auf, 0 
läßt sich leicht die dritte Dezimale von f ablesen, - 
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werte und Lösungen der Integrodifferentialgleichung 
(1) ziehen kann, sei hier verzichtet. 

Für die praktische Berechnung von Eigenwerten 
der Integrodifferentialgleichung (1) ist von Nutzen, 
daß man die für Integralgleichungen entwickelten 
Näherungsverfahren (wie z.B. das Schwerinsche oder 
Löschsche Verfahren)!) nunmehr auf die begleitende 
Integralgleichung anwenden kann. 


Stuttgart-Vaihingen. W. Buscham. 


1) Siehe Z. angew. Math. Mech. Bd.24, 8.35—41 und 
2. techn. Phys. Bd.8, 8. 264—271. 


Die Potenzreihen der Ogivalfunktionen. 
„‚Ogival‘“ nennt man den Drehkörper eines Spitz- 
bogens. Ein Ogival entsteht also durch die Um- 
drehung eines Kreisbogens vom Radius R und Win- 
kela um eine Achse, welche durch einen Endpunkt des 


Kreisbogens geht und der Tangente im anderen End- 
punkt parallel ist. Das Doppelte des größten Abstan- 
des des Kreisbogens. von der Achse gibt dann den 
Durchmesser D des Ogivals, genauer: seines Grund- 


kreises. Die Spitze des Ogivals befindet sich in einer 
Höhe H über dem Grundkreis; der halbe Winkel an 
der Spitzeist gleich a. R heißt der Abrundungsradius 


des Ogivals, das Verhältnis > =n das Abrundungs- 


verhältnis oder die Abrundungszahl. n bestimmt die 
Form des Ogivals, D (oder auch R oder H) die abso- 


lute Größe. Der Falln = 4 entspricht der Halbkugel. 


In den Formelsammlungen, welche Volumen, Ober- 
fläche, Schwerpunkt, Trägheitsmomente geometrisch 


einfacher Körper enthalten, finden sich manchmal 


(allerdings nur selten) auch Angaben für das Volumen 
V und die Mantelfläche M eines Ogivals. Doch sind 
die angegebenen Formeln, außer für kleinere n (bis 
etwa n = 3), für die Auswertung unbequem. Wohl 


führen alle die hier in Betracht kommenden Integrale 
nur auf elementare Funktionen. Aber man erhält z.B. 


für das Volumen die Formel 


1% 


"und ist nun etwa n = 6,5, wobei H = 2,5 D wird, so 


kommt 


Er Y=,rD® (1205 — 30420). 


Da nun in diesem Beispiele die Klammer den Wert 
4,0456 hat, so muß man & mit relativ großer Genauig- 
keit bestimmen; der dreistellige Wert x = 0,395 liefert 

nur 3,410 statt 4,0456. j j 
Esist daher vorteilhafter, falls man öfters mit Ogi- 
valen zu tun hat, den asymptotischen Wert 
si . % D} h O% . 


N>& 


heranzuziehen und einen Korrekturfaktor f(n) hinzu- 


die gut konvergente Reihenentwicklung nach fallen- 


mit V hinreichend genau zu berechnen ist. — Für die 


noch in verstärktem Maße zu. 


Setzt man H= Dh, so hängt h nur von n ab: 


a 


Br: A sin = + 
0088-15, Be 2 2/n 


D2[H(12n®— An +1) —-Dal2m®—6n)], 


zufügen, der für n— oo gegen 1 strebt und sich durch 
Momente von V und M treffen solche Erwägungen > 


Winkel « bestehen die leicht ersichtlichen Be- 
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Be ee Ve a ET Tg 


Damit wird die oben angegebene Volumenformel im 
allgemeinen Falle 


2 __%_(19n2 6 | 
vV= men | 12n® —4n +1) se (12% n) 
— Lupe zn? ni ya a2n3)|. 
12 tga 


Setzt man hier die Reihen ein, die sich auf Grund der 
Beziehung zwischen & und n ergeben: 


Dre ee a 
en = pr 
SIE: ; 
7 an, 
Fe) (2i+1)! 
& ET RRTD 3 
en dB 
SG le +)! 
—=-1— VEIT IIE NE -—i 
2 it! "> 


so entsteht die Entwicklung 


ri 2 SDü+dldtg! |, 
v-amı (+3 @ir+H! =) 


2 
— 757 D’hfm, 


woraus die Reihe für die Ogivalfunktion f(n) unmittel- 
bar abzulesen ist. In entsprechender Weise sind die 
Reihenentwicklungen .für die 
anderenOgivalfunktionen g (n), 
k(n), Un), m(n), p(n) gewonnen 
worden. In ein Kurvenblatt, 
der Ogivalfunktionen nimmt 
man zweckmäßigerweise auch , 
& und h(n) auf, für welche 
keine Reihenentwicklungen in 
Frage kommen. 

In den nachstehenden For- 
meln bezeichnet V das Vo- 
lumen, M die Mantelfläche Schwerpunktshöhe 
eines (durch D und n gegebe- 
nen) Ogivals; unter Voraus- 
setzung einer gleichmäßigen 
Raumerfüllung Fer Flächen- 
belegung mit Masse sind 
dann jeweils (d.h. für V bzw. 
für M) x die Höhe des Schwer- 
punktes über der Ebene des 
Grundkreises, r der Trägheits- 
radius für Drehungen um die 
Ogivalachse (der ‚‚Drallarm‘“), 
ge der Trägheitsradius für 
Drehungen um einen Grund- 
kreisdurchmesser (der ‚‚Tau- 
melarm“). Wird das Trägheits- 
moment für Drehbewegungen um eine Achse (also: die 
Drehmasse um diese Achse) gesucht, die im Abstande a 
unterhalb der Grundkreisebene parallel zu dieser durch 
die Verlängerung der Ogivalachse geht, so ist 0? zu 
ersetzen durch (9? + 2ax + a®). Insbesondere gibt 
@ = — x das äquatoriale Hauptträgheitsmoment (die 
Haupttaumelmasse) des Ogivals. Doch da in den An- 
wendungen das Ogivalin der Regel nur einen Bestand- 
teil einer zusammengesetzten Körperform ausmacht, 
so ist dieses auf den Schwe akt des Ogivals be- 
zogene Moment hier nicht aufgeführt, 


Reihen der Ogivalfunktionen 
Volumenfaktor f: 


Größe 


* Schweremoment 


Drallmasse 


(Drallarm)? 


Taumelmasse 


(Taumelarm)? 


1 


1 
De Aggntt nt IE a 


336 308 6864 


an) = — Ir Br BR +19: + 24) ni 
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Mantelfaktor g: 
3 ee A ee 
g=l +59” +35” + og” + 70” 


Drallmassenfaktor des Volumens k: 


1 — 21 -2 29% 83 543 _ 
kalt FT ao aan 


Taumelmassenfaktor des Volumens I: 


1 19 ". 203 1 
per 1 ___m-2 anne en ... 
i=1l +7n +95” +75" +94” + 


Drallmassenfaktor des Mantels m: 
5 41 49 1 
2 Dr ER RE 
Tas OuB2 2 ae 
Taumelmassenfaktor des Mantels p: 
1 2a u 2008 . 83 
EINE Re RE BAER. 
Pl Far Tao ara 
Zum Schluß mögen noch die Ausdrücke für die all- 
gemeinen Reihenglieder angeschrieben werden: 


15 (i+ 1 +3)! 
ar = 


Aloe. 


n% + ... 


FIT BR 
on) Haren! De 


16 @i+9 


kin) =1 + DJ rss) 
i=1 


Formeln für das Ogival 


6 m 
SIE ech 
r 385 D? - 
Mo? == = rn Dinh 
15 2 
" 1 np 
8.0 —. Da. 
J N 5” 9 
(x und o sind auf einen Grundkreisdurchmesser bezogen). 
In)’ =1+ 2 I 


4 
I 105 id +4)! Y ’ 
2 16 @irrı (9E+ 358 4 46)m 


135 il(i +)! 


Ki; I 15 (di Didi+3)! 
p(n) = A 2 1 its (742 + 11i+10)n-i. 
Die Reihen konvergieren absolut für n > n En 


Brunoy. U.T.Bödewadt. 
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Bemerkungen zur numerischen Lösung von 
Randwertaufgaben für nichtlineare gewöhn- 
liche Differentialgleichungen nach der Picard- 
schen Iterationsmethode. 


I. Einleitung. 


In einer früheren Arbeit!) wurde eine Methode zur 
numerischen Lösung von Randwertaufgaben für ge- 
wöhnliche lineare Differentialgleichungen entwickelt. 
Diese Methode, die mit Reihenentwicklungen nach 
Legendreschen Polynomen sowohl für die Koeffizien- 
ten als auch für die gesuchte Lösung der linearen Diffe- 
rentialgleichung operierte, gestattete die Berechnung 
der allgemeinen Lösung der vorgegebenen Differential- 
gleichung. Dies war ein gewisser Vorzug der Methode, 
der z.B. dann zur Geltung kommt, wenn die Rand- 
werte nicht von vornherein festgelegt werden, sondern 
noch offen bleiben sollen. Unter diesem Gesichtspunkt 
erschien der Rechenaufwand dieser Methode (Auf- 
lösung von sechs linearen Gleichungen) tragbar. 

Wenn man indessen auf die Gewinnung der all- 
gemeinen Lösung der Differentialgleichung verzichten 
kann und sich auf die Behandlung einer Randwert- 
aufgabe mit fest vorgegebenen Randwerten be- 
schränkt, wie man dies auch bei den bekannten Lö- 


- sungsverfahren im allgemeinen tut, kann man diese 


spezielle Lösung unter Verwendung der Legendre- 
schen Entwicklungen nach der Picardschen Iterations- 
methode mit im allgemeinen ziemlich geringem 
Rechenaufwand gewinnen. Die Verwendung der 
Legendreschen Entwicklungen gestattet hierbei eine 
bequeme Berechnung der einzelnen Iterationen. 

Dies soll im folgenden an Hand der ersten Rand- 
wertaufgabe für die Differentialgleichung zweiter Ord- 
nung gezeigt werden. Die Beschränkung auf lineare 
Differentialgleichungen, die in der oben zitierten Ar- 
beit noch notwendig war, ist hier natürlich nicht er- 
forderlich. 


II. Die Picardsche Iterationsmethode unter 
Verwendung von Reihenentwicklungen nach 
Legendreschen Polynomen. 

Die Differentialgleichung zweiter Ordnung sei vor” 
gegeben in der Form: 


SERIEN SB er (1). 


Die Funktion 9 möge allen Bedingungen genügen, die 
erforderlich sind, um die eindeutige Lösung der ersten 
Randwertaufgabe: 


ee 
z«=—1l:y=y(—]) 


sicherzustellen?). . 

Bei der Picardschen Iterationsmethode wird in die 
rechte Seite von Gl.(1) eine Näherung der Lösung 
eingesetzt, die den Randbedingungen Gl.(2) genügt, 
und aus Gl.(1) dann eine verbesserte Näherung ge- 
wonnen, Durch das Einsetzen der Näherung in die 


- rechte Seite von Gl. (1) wird diese eine Funktion allein 


von %&: 

EN Neo ae EEE (3). 
Die rechte Seite von Gl.(3) entwickle man jetzt nach 
Legendreschen Polynomen Py(z)?): 


Je) = 2, Pole) rs (E) 


“und mache auch für y(x) einen solchen Entwicklungs- 


ansatz mit unbestimmten Koeffizienten Y»: 
; ya)=&XYor- Pole). - -- (5). 
v ’ 


7Z. w.Math. Mech. 31 (1951), 8. 104/114. 
s Denen die Randwerte an den Stellen © = a und & = b, so 
bringe man diese durch die lineare Transformation: E 
ed. 5,222 ee. 


3) Solche Bedingungen sind angegeben z. B. bi Kamke: 


" Differentialgleichungen, Lösungsmethoden und Lösungen, Bd.1, 
ATeipaig 1944, 8. 277 ff 
) Zu 


r Berechnung der Entwicklungskosffizienten vgl. Z, angew. 


Math. Mech. 44 (1944), 8. 71/76, 


Setzt man Gl.(4) und (5) in Gl.(3) ein und drückt 
Pj(z) wieder durch die Legendreschen Polynome 
selbst aus, so kann man aus Gl. (3) durch Koeffizien- 
tenvergleich der Legendreschen Polynome ein lineares 
Gleichungssystem für die gesuchten Koeffizienten Yy 
gewinnen, das als rechte Seiten die (bekannten) Ent- 
wicklungskoeffizienten F, von ffz) enthält. 

Für den K.oeffizientenvergleich denkt man sich die 
Reihen Gl.(4) und (5) abgebrochen, sobald die Ent- 
wicklungskoeffizienten hinreichend klein geworden 
sind. Dies ist im allgemeinen der Fall, wenn man » 
maximal bis 7 laufen läßt. In diesem Fall erhält man 
aus Gl.(3) durch Koeffizientenvergleich von P,(z), 
Pı(z), Ps(&), Pz(x), Pı(x) und P,(x) das folgende Glei- 
chungssystem: 


sr. Eur Formen 
15 Y, +42 9,7 sı.9, =7] 
Br + WwH=r, 


ER De (6) 
9Y,=F, 
143 Y,=F, 


Zu diesen sechs Gleichungen für die acht Koeffizienten 
Yo I‘, Y, treten jetzt noch aus den Randbedin- 
gungen Gl.(2) die beiden Beziehungen: 
Y+MH+B+ YAM YrtYtH=yH) (7). 
Y .utrh, Mile hut Henn) 


Mit den Abkürzungen: 


Zt +n]=a = 
+ snI=R | 


erhält man aus den beiden Gl. (7) durch Addition bzw. 
Subtraktion: 


ren. 
+, +, rF9=B 


Das aus Gl. (6) und (9) zusammengesetzte Gleichungs- 
system läßt sich, von Y', bzw. Y, ausgehend, sukzes- 
sive auflösen. Als Lösung ergibt sich: 


7 1 1 
Y=&- Ft 7zF: 


ee ; 
D=,% antun 
er ar ... 0) 
Kanon 
Em r 

K=a4r 

= 


Damitist die nächste Näherung Gl. (5) des Iterations- 
verfahrens gewonnen®). Diese hat man wieder in die 
rechte Seite von Gl. (1) einzuführen, die dadurch ent- 
stehende neue Funktion f(x) wieder nach Legendre- 


>) Zur Gewinnung der ersten Näherung setze man z.B. alle 
F,in @1,(10) gleich Null. Dann erhält man als erste Näherung 
die geradlinige Verbindungsstrecke der beiden Randpunkte (—1, 
y(-2) und (+1,y(+ 1)). 
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schen Polynomen zu entwickeln und mit diesen neuen 
Koeffizienten F, erneut Gl. (10) zu bilden usw.°)?). 

Bei Differentiagleichungen von höherer als zweiter 
Ordnung läßt sich die Picardsche Iterationsmethode 
in ähnlicher Weise durch Legendresche Entwicklungen 
numerisch bequem durchführen. Da hierfür die Ab- 
leitungen höherer Ordnung der Legendreschen Poly- 
nome, durch letztere selbst ausgedrückt, benötigt wer- 
den, seien diese Ableitungen bis zur vierten Ordnung 
hier noch angegeben ®): 


P,)=(2n—1)  Pa-ı + (2n—5) Ps 
+ (2n—9) Pn-5 + (2n—13)  Pn-„ (n 7) 
P} = (2n—1)(2n—3) P)-3 + 2(2n—3) (2n—7) 
Pn-ı + 3 (2n—5) (2n—11) Pn-; (nSS7) 
P% = (2n—1) (2n—3) (2n—5) Pn-3 
+ 3 (2n—3) (2n— 5) (2n—9) Pan-5 
+6(2n—5) (2n—7) (2n—13) Pn-, (n<7) 
Pr"—(2n—1) (2n—3) (2n—5) (2n—7) Pn-ı 
+4(2n—3) (2n—5) (2n—7) (2n—11)Pn-s 
(n ST) 


(11). 


III. Zwei Beispiele 


Diese hat als Lösung bekanntlich die Kettenlinie: 
y= 16022 — &0j2} ru 


Beginnt man das Iterationsverfahren mit der gerad- 
linigen Verbindungsstrecke der Randpunkte, d.h. in 
diesem Falle mit der Strecke -1<rs-+l1 der 
x-Achse, so erhält man für die Legendreschen Koeffi- 
zienten Yy der nächsten Näherungen die folgenden 
Werte der folgenden oberen Tabelle. 

Die folgende untere Tabelle bietet für einige x-Stel- 


len einen Vergleich der Näherungen y;=% 7% -Py(z) 
v 


mit der exakten Lösung ynach Gl. (13): 
An den ausgerechneten Stellen weicht also y, von 


y um höchstens 3 oo ab. 


In diesem Beispiel hat man nach jedem Iterations- 
schritt die rechte Seite f(x) von Gl.(3) numerisch 
nach Legendreschen Polynomen zu entwickeln. Dies 
ist indessen nicht immer erforderlich, vielmehr kann 
man bei manchen Differentialgleichungen wenigstens 
die ersten Iterationsschritte noch einfacher gewinnen. 
Man vergleiche hierzu das folgende Beispiel: 

b) Es sei die erste Randwertaufgabe: 


Amar 1 
a) Es sei die erste Randwertaufgabe: 2 = —1, y=0, i=—1l, y=--; ı=+l y=— 
z=-+1, y=0 für die Differentialgleichung: 2 4 
y'=2: yı+y® RT <T (12) für die een ; 2 
. = — yy N ae ae 
gestellt. gestellt. 
I. i=2 i=3) | iu4 i=5 i=6 i=7 
r® Fa} —0,87261 —0,94654 —0,96801 —0,97310 —0,97410 
r 2, 0,82041 0,86579 0,87685 0,87909 0,87950 
ri 0 0,05450 0,08057 0,08885 0,09072 0,09099 
rW 0 —0,00230 0,00018 0,00231 0,00329 0,00361 
RE WER ABLE SBERSEEBRRERREN 7 SER BEE ENT IE IT 
c | Yy Ya Ya Yı Vs Ys Y%ı 
() —1,3811 —1,0000 —1,2617 —1,3493 — 1,3738 -—1,3797 —1,3809 
0,2 —1,3406 —0,9600 —1,2208 —1,3088 — 1,3334 —1,3391 — 1,3403 
0,4 —1,2124 — 0,8400 —1,0927 —1,1807 — 1,2054 —1,2110 —1,2120 
0,6 —0,9758 — 0,6400 —0,8624 | —0,9447 —0,9688 — 0,9744 — 0,9754 
0,8 —0,5924 —0,3600 —0,5070 —0,5671 —0,5863 —0,5911 —0,5921 
1,0 0 0 0 0 0 0 0 


®) Liegt nicht die erste, sondern eine andere Randwertaufgabe 
vor, so ändern sich in G1.(10) nur die Ausdrücke für die Koeffi- 
zienten Y, und Y,. Auch die Anfangswertaufgabe kann man, falls 
man die Lösung nur für ein endliches «-Intervall benötigt, ganz 
analog behandeln. Dieses endliche x-Intervall bringe man durch 
eine lineare Transformation auf das Intervall -1IsSa<s +1, 
Für die Anfangswertaufgabe: x = — 1, y = y(—1), y’ = y’(—1) 
sind die beiden ersten Gleichungen von (10) zu ersetzen durch: 


n y 2 
a Tu VE us u 2 ee ee Ze 


Hey(-l)+ n—4n+ ER 


während die übrigen wieder ungeändert bleiben. 


”) ImSpezialfall der linearen Differentialgleichung: y’’+a(z)y’ 
+b(2)y=c(e)kannmannach Picard das Iterationsverfahren 
auch so durchführen, daß man zuerst die Randwertaufgabe Gl. (2) 
für die Differentialgleichung: y’’= c(«) löst. Die durch G1. (10) 
gegebene Lösung, bei der jetzt F, die Legendreschen Koeffizienten 
von c(®) bedeuten, möge y, genannt werden, Als nächsten Schritt 
hat man dann; y’’= a (2)y, + biz)y für die Randwerte y(+1) 
=y(—1)=0 zu lösen. Die Lösung sei y,. Diese setze man 
; dann statt y, in die letzte Differentialgleichung ein, bestimme 
daraus wieder für die Randwerte y(+1) =y(—1) =0 die Lösung 
Ys nach G1.(10) usw. Die gesuchte Lösung y(x) erhält man dann 
aus Y(@) = Yol®) — Yılz) + Yalz) — Yalz) + —'''. Diese für 
lineare Differentialgleichungen anwendbare Form der Iteration ist 
numerisch besonders angenehm zu handhaben, da die Koeffizien- 


ten ri) von y;(z) mit wachsendem ö immer kleiner werden. 
®) Glieder mit P,(z),in denen y negativ ist, sind wegzulassen. 


” 


Diese Aufgabe hat, wie man leicht nachprüft, die 
Lösung: | 


re 

fear oe rs 

Die Legendreschen Koeffizienten F, der rechten 

Seite von Gl.(14) nach Einsetzen einer Näherungs- 

lösung y(z)=% Yy Py(x) kann man z. B. für die ersten 
v 


drei Iterationsschritte gewinnen, indem man in der 
letzten Summe » von O0 bis 3 laufen läßt. Drückt man 
nämlich die Produkte P,(z) P’ (x), die nach dem 
Einsetzen auf der rechten Seite von Gl.(14) ent- 
stehen, wieder durch die Legendreschen Polynome 
selbst aus, so findet man für die F, im Falle der Dif- 
ferentialgleichung (14): 
R=-2-(YeH+Hh+YYı+Y,%) 

F. 1 2 


(Ve HrHBrer+ ger 


(16) 
F,=—2: (rer nr) 


n=-2:- Bazz: tz) 
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Er mas, ur, | mit der Funktion f(x) übereinstimmt, und auch die 
® er, gesuchte Iteration als Parabel vierter Ordnung ansetzt. 
2: 4 or (16).. Da dann Y,=Y,=Y,=0 ist, fallen jetzt die drei 
FR=-—r:. letzten G1.(6) fort. An Stelle von Gl.(10) erhält man 
21 p jetzt die folgende Auflösung: 
en: 2 : 3 1 z 1 1 
Die Ausdrücke Gl1.(8) werden hier: «= 3 ß ne Y,=a— #0 +7zFr 
Aus Gl. (10) findet man dann für die erste Näherung, 1 
indem man alle Fy, gleich Null setzt, Y, =ß-— is? 1 
3 1 
rw=-, rw=_—, Y=- Yb=-.=r70- 
en z. Yü=Y0 Yvv=0. Ges sr, il: 
Mit diesen Werten folgt aus Gl.(16): #,=0,093750, 
F,=—0,03120, P,=F,=F,—F,-0 und aus ven 
Gl.(10) für die zweite Näherung: 15 
Y@ = 0,34375,  Y{)= —0,12292, a n F, 


Y@ = 0,03125, Y@ = —.0,00208. 

Diese Werte kann man nun nochmals in Gl. (16) 
einsetzen, da in G1l.(16) gerade die Entwicklungs- 
koeffizienten bis Y, einschließlich berücksichtigt wor- 

_  densind. Man erhält so neue Werte für FF, **', Fy. 

Aus Gl.(10) findet man dann daraus als dritte Nähe- 
tung: ; : 

79 =0,34580, Y@ =—0,11859, Y® =0,02834, 


Y@) =-+0,00001, 7 =0,00000. 


_ Die folgende Tabelle zeigt für einige x-Werte die 
zweite und dritte Näherung im Vergleich zur exakten 
Lösung y: Ba 


z IB Bde | % %» 
Er 0,5000 


Be ic 0,5000 | 0,5000 

0,6 | 0,4167 | 0,4180 | 0,4154 

0,2 | 0,3571 | 0,3540 | 0,3555 

0,2 | 0,3125 | 0,3060 | 0,3116 
0,6 | 0,2778 | 0,2720 | 0,9777 

1 - | 0,2500 | 0,2500 | 0,2500 


t diese Genauigkeit noch nicht ausreichend, so 
müßte man für die weiteren Näherungen jetzt jedes- 
mal die rechte Seite f(x) von Gl.(3) numerisch nach 

chen Polynomen entwickeln. 


tberücksiehtigen. Diese werden indes- 
dem » immer umfangreicher, so daß 
1. (14) nach jedem Iterations- 


ıtwicklungsverfahren 
nfür dieersten 


in III. legt 


en zu suchen, 
lgleichungen 
auf jede 


7, 


.  7@=--0,00632, Y@) —=0,00086, Y7@—=—0,00009, 


ı könnte natürlich zu Gl. (16) analoge Ausdrücke 
ie Fy aufstellen, die auch die Glieder Yy für 


r gegenüber der Neuentwicklung - 
G 


tionsschritte 


einen E 


Es’bleibt noch die Berechnung der Koeffizienten F„ 
für die ersten Iterationsschritte zu beschreiben. Diese 
erfolgt nach G1.(5) der oben zitierten Arbeit!), in der 


jetzt | Ar] = 7 zu setzen ist. Für F„ ergibt sich hier- 


aus dann: 


<0;17 
1 : 


me, + Dr en +Dnon | 


+4: [Afe HD" Apal- On + Ku] 
+8- [fe HD*- A] (2m +1) 700,1) 
+ AH Hr 20 | 
E - (@n+1) Mn(0,1) 
+ [Atfa+ (Ir Ara] 
| -OR+DN.OD) 


An die Stelle der Tabellen 1—5 der oben zitierten 


Arbeit!) tritt jetzt die folgende Tabelle, die man aus 
- den Definitionsgleichungen für die In, Kn; In, Mn» 


N,„ direkt berechnen kann: 


700,1) | Kn(0,1) | Ln(0,1) | Mn(0,1) | No) 3 
Be > ae 
£ 2]sr 2 24 16 | 960 
je N 
REN 2 50 320° 
En TR BR 181 
a Re By 56 a | 3% 


 — Nach Bildung der Funktionswerte f(x) anden Stellen 


en NO RIE wie 
ı=0, Zr Eior, > | 


ten Differenzen einschließlich) hat man dann (nach 


der Faltung) lediglich 5 Multiplikationen unter B- 
nutzung der obigen Tabelle auszuführen, um -— in 
einem Rechengang auf der Rechenmaschine — je 
ntwicklungskoeffizienten F„ zu erhalten. Bei 

den oben zitierten Entwicklungsverfahren )!) 


erfür 15 bzw. 25 Multiplikationen erforderlich. 
ptsächliche Arbeitsersparnis liegt jedoch bei 


n Differentialgleichungen (1) kann die 

Funktionswerte ziemlich zeitraubend 
m soeben beschriebenen vereinfachten 
sverfahren für die ersten Iterations- 
diese Funktionswerte nun nur an 9 Stellen 
ntervalls zu bilden gegenüber 21,bzw. 41 Stellen 


Re . e F 


Je 


£ gr f DEV: 
a > er 


ne On+1).|On +1). 


und Aufstellung des Differenzenschemas (bis zu vier- 


der Funktionswerte f(x). Bei etwas 
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bei Zugrundelegung des genaueren Entwicklungsver- 
fahrens unter Benutzung von Näherungsparabeln 2ter 
bzw. 4ter Ordnung. 

Für die höheren Iterationen wird dieses vereinfachte 
Entwicklungsverfahren indessen im allgemeinen nicht 
genau genug sein, da es für die Iterationen Gl. (5) 
immer nur Parabeln 4ter Ordnung liefert. 

Die Anwendung dieses vereinfachten Entwicklungs- 
verfahrens auf das Beispiel IIIa ergibt für die Ent- 
wicklungskoeffizienten der ersten Iterationen die 
folgenden Werte: 


yo |— 3 —0,87247 | —0,94929 | —0,96995 
y®l+ 3 +0,82387 | -+0,86909 | +0,87497 
y@| 0| -+0,04859 | -+0,08020 | +0,09498 
Frankfurt/Main. E. Fehlberg. 


Einfache Herleitung der allgemeinen 
Kirchhoffschen Beugungsformel und ihres 
elektromagnetischen Analogons. 


Die vorzügliche Monographie über das Huygensche 
Prinzipvon Baker und Copsonist mit dem 
Erscheinen der zweiten Auflage auch in Deutschland 
zugänglich geworden. Bei Durchsicht des Buches fällt 
auf, daß — im Gegensatz zu der eleganten Durch- 
führung vieler anderer Ableitungen — die allgemeine 
Kirchhoffsche Beugungsformel sehr umständlich 
und unübersichtlichhergeleitet wird,unter Aufwand von 
über zwei Seiten mit mehr als zwanzig Formelzeilen!). 
Es scheint demnach weder allgemein bekannt noch 
trivial zu sein, daß diese, Gleichung in wenigen Zeilen 
gewonnen werden kann, nämlich im wesentlichen 
durch nichts weiter als eine zweifache partielle Inte- 
gration. Diese einfache Ableitung sei im folgenden 
mitgeteilt (implizit ist sie bereits in einer früheren 
Arbeit des Verfassers enthalten?). Darüber hinaus 
wird anschließend die analoge Ableitung für elektro- 
magnetische Wellen und für höhere Räume ungerader 
Dimension gegeben. 

kustische Wellen. Statt von der homo- 
genen Wellengleichung gehen wir allgemeiner von der 
inhomogenen aus: 


1.9 
Az =—0 Re sei (1). 


Ohne daß der Gang der Rechnung durch diese Ver- 
allgemeinerung geändert wird, erhalten wir so mit 
einem Schlag die Kirchhoffsche Beugungsformel wie 
auch die retardierte Lösung der inhomogenen Wellen- 
ns: — Wir betrachten ein Volumen V des 
Raumes, nennen t den Radiusvektor eines Punktes J 
in V, und r den Abstand des Punktes J von einem 
Aufpunkt P außerhalb VY. Unter dem retardierten 
Wert einer Funktion F(r,t) des Radiusvektors und 
der Zeit verstehen wir 


F@I=rln-n): ur Aa 


Für eine solche Funktion gilt ersichtlich nach dem 
verallgemeinerten Greenschen Integralsatz 
angewandt auf die Funktion x(r) [F] 


S xt) [WFJde = —f do xe)[F]+ VS Kr) LF]de (3). 


ı)B.B. Baker und E. T. Copson, The Mathematical ° 


Theory of Huygens’ Prineiple. Oxford University Press, 2. Aufl. 
1950 (1. Aufl. 1939), Kap.I, $ 5.2. 

®) W. Franz, On the Theory of Diffraction. Proc, Phys. 
8oc. London A 63, 8. 925 (1950). 


Z. angew. Math, Mech. 
Bd. 32 Nr.1ı Jan. 1952 


Darin ist dr das Volumelement, do das vektorielle 
Oberflächenelement von V (do habe die Richtung der 
inneren Normalen), x eine differenzierbare Funk- 
tion von r; die Funktion F kann wahlweise ein Skalar 
sein oder ein Vektor, welcher skalar oder vektoriell 
mit y bzw. do multipliziert ist. Die Integrationen 
sind über das gesamte Volumen V bzw. dessen Ober- 
fläche zu erstrecken. Ferner wollen wir verabreden 
(um Indizes zu vermeiden), daß Differentiationen unter 
dem Integral stets nach dem Radiusvektor des Inte- 
grationspunktes J ausgeführt, solche nach dem Auf- 
punkt P dagegen vor das Integral geschrieben wer- 
den. — Wir setzen y=1/rund YF=Y grad u, und 
erhalten durch nochmalige Anwendung von Gl.(3) 
mit f=u: : 


dt ı1 024 do 
[Fe-s#|-- [Fee 
air (Frt (az) | 


Der letzte Summand verschwindet, da sein Inte- 
grand die Gestalt f((—r/c)/r hat und daher der homo- 
genen Wellengleichung genügt. Ist u Lösung der 
Wellengleichung (1), dann können wir die linke Seite 
mit ihrer Hilfe umrechnen und erhalten 


[F@- Grad] dir | 7-0 (5). 


. (4). 


Dies ist bereits die Kirchhoffsche Formel (vermehrt 
um den ersten Term) für den Fall, daß der Aufpunkt P 
außerhalb V liegt, wie wir bis jetzt angenommen 
haben. Liegt P innerhalb V, dann schließt man ihn 
in üblicher Weise durch eine kleine Kugel aus, deren 
Radius man gegen Null gehen läßt. P bleibt auf 
diese Weise außerhalb des Integrationsgebietes, so 
daß G1.(5) gilt; jedoch liefert die Oberfläche der un- 
endlich kleinen Kugel einen Beitrag zum dritten Sum- 
manden von Gl.(5), und zwar, wie man unmittelbar 
sieht, —4ru(P). So ergibt sich für P innerhalb V: 


d 
inı= (Ftl- _ Ierad u] - iv [Lu (6). 


Diese Gleichung vereinigt die Kirchhoffsche Beu- 
gungsgleichung (für e=0) mit der Ausstrahlungs- 
Lösung der inhomogenen Wellengleichung für den un- 
begrenzten Raum (hierbei verschwinden die Ober- 
flächenintegrale). 


Elektromagnetische Wellen. Wir gehen 
aus von den Maxwell-Gleichungen 


1069 =e%; dv9=0; 
cl 
rt En; dive=0 |. 


Durch zweimalige partielle Integration ergibt sich 
gemäß Gl.(3) (mit VF=P xrot €, beim zweiten Mal 
FV=Vx6) 


dr ı € do 
j 2 [rot rot C+ = 3 =— 7 x [rot €] @ 
do 1.92 dt 2 
— rot j „xee + (rotrot = Fa) j > [E] 


Die linke Seite verschwindet nach Gl.(7), und der 
letzte Integrand rechts genügt als f(t—r/c)/r derhomo- 
genen Wellengleichung, und daher ist wegenrotrot - 
= graddiv —4: 


(rotzot + zn) | 2 [6] = grad div f = [6] 


(9). 
= grad j 7 De] 


E worin die Funktion © des Differentialoperators 7 12 
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Der letzte Schritt folgt mittels G1.(3) aus divE=0. 
Aus Gl.(8) ergibt sich somit 


Frxtorei+ror = x [€] | 


(10). 
ana | F-[E1=0 


Verlegen wir den Aufpunkt P ins Innere von V und 
schließen ihn durch eine unendlich kleine Kugel aus, 
dannliefert diese nur zu dem zweiten Summanden von 
Gl.(10) einen endlichen Beitrag (da nur in diesem 
Integral ein gerader Anteil —1/r? enthalten ist), und 
zwar —4nC(P), und es ergibt sich die Gleichung 
vonLarmorund Tedone(s. Baker-Cop- 
son, 8.102ff.) 


a EST SLIEE x[€] 
ana | 7-16) 


Diese Gleichung bestimmt den Wert von E in F 
aus den Randwerten (die analoge Gleichung gilt natür- 
lich auch für 9), leistet also für die Maxwellgleichungen 
dasselbe, wie die Kirchhoffsche GI. (6) für die Wellen- 
gleichung. Dennoch besteht für die beugungstheo- 
retische Anwendung ein wesentlicher Unterschied 
zwischen G1.(6) und (11): G1.(6) liefert für beliebige 
(also auch angenäherte) Randwerte Lösungen 
der akustischen Wellengleichung, da die Integranden 
diese einzeln erfüllen. Die Wellengleichung wird zwar 
auch von Gl.(11) erfüllt, jedoch nicht die Bedingung 
div&E=0, welche für Lösungen der Maxwellschen 
Gleichungen zusätzlich zu fordern ist. Daher ist 
Gl.(11), im Gegensatz zu Gl.(6), für Zwecke der 
Beugungstheorie ungeeignet. Man kann jedoch diesen 
Übelstand beseitigen und die Divergenzbedingung 
identisch (d.h. für beliebige Randwerte) befriedigen, 
indem man von Gl.(11) die Rotation bildet. Unter 
Benutzung von Gl.(7) entsteht dann 


(11). 


d d 
ee) _ ek EG (12). 


Die entsprechende Gleichung für die elektrische 
Feldstärke ee 


Anz 0 


1 d 
En a [7 [€] + rotrot 7% [9] (12a) 


_ erhält man, indem man in G1.(11) & durch $ ersetzt 


und dann die Rotation bildet, oder auch, indem man 
diese Gleichung nach der Zeit differenziert (dabei 
bringt man das letzte Glied durch eine zu Gl.(3) 
analoge Anwendung des Stokesschen Integralsatzes 
zum Verschwinden). — Die Gl.(12) und (12a) liefern 
nicht die Feldstärken selbst, sondern .deren zeitliche 
Ableitungen; für diese gelten jedoch die Maxwellschen 
"Gleichungen identisch, 


Akustische Wellen im m-dimensionalen 
Raum (m ungerade > 3). Die Rechnung für den 


- dreidimensionalen Raum überträgt sich fast ungeän- 


dert auf höhere Räume ungerader Dimension. Der 
"einzige Unterschied besteht darin, daß nicht f(1—r/e)/r 
der homogenen Wellengleichung genügt, sondern statt 
dessen gilt (wie man unmittelbar verifiziert) 


Ba rd r 
(4 Sue 2(a)t | = See 


co 
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das folgende Polynom ist: 


Se 2 
v 


D(e) = 2, (22)” er 
> ine” ® (&)? "mu 2 
m—3 


Der Polynom-Ansatz Bist nur möglich, wenn 


eine ganze Zahl, also die Dimension m ungerade ist. 
Deshalb beschränken wir uns auf diesen Fall, ohne 
auf die bei gerader Dimension eintretenden Verhält- 
nisse einzugehen?), — Wir können jetzt die Gl.(4) 
und (5) übernehmen, wenn wir darin nur unter den 


Integralen 1/r durch = 6) | 0 ersetzen. Wird der 


\ 
cä ) 
Aufpunkt P ins Innere von V verlegt, dann liefert die 
P ausschließende Kugel wieder nur einen Beitrag zu 
dem dritten Summanden von Gl.(5), welcher aber 
jetzt ersichtlich gleich — (m — 2)0mu wird, worin om 
die Fläche der m-dimensionalen Einheitskugel ist. 
An Stelle von Gl.(6) ergibt sich also 


m don = | 0) 


-® © =) [grad u] ! (15). 


: do ) 
= div j rm—2 ® FE a) [%] 


Walter Franz. 


Münster. 


a) S. Courant-Hilbert, Methoden der mathematischen 
Physik, Bd. 2 (1937), 8. 395. 


Über iterierte Integration von Bessel- 
Funktionen. , 
(Mitteilung aus dem Heinrich-Hertz-Institut) 
@ 


1. Einleitung 

Bei verschiedenen Problemen der Theorie von Ein- 
schwingvorgängen in Kettenleitern spielen mehrfache 
Integrale über Bessel-Funktionen eine Rolle. 

Es zeigtsich, daß auch hiersich Rekursions-Formeln 
aufstellen lassen, in denen als wesentliche Funktionen 
außer den Funktionen Nullter Ordnung J,(x) und 
erster Ordnung J, (x) das erste Integral über J,(2) 


% 
SIolt)dt 
ö 
auftritt. 
© 
Das Integral /Jo(t)dt ist weit weniger numerisch 
0 


erfaßt als die Besselschen Funktionen J, (x) und J; (2). 
Es befindet sich tabelliert in Watsons Werk ‚,‚A trea- 
tise on the theory of Bessel-Functions‘; jedoch nur 


für Argumente x zwischen 0 und 1. 


Zur numerischen Berechnung kann außer Reihen- 
Entwicklung verwendet werden die Beziehung (vgl. 
Watson, Kap. 16,31): 


F In) —2r SZ Iptsatı(e) ERDE 
0 n=0 


Ferner _die Beziehung, welche auf die Struvesche 
Funktion H,(®) und ihre erste Ableitung zurückgeht 
(vgl. Watson, Kap. 10, 74): 


[ Jet = Wola)- Hua) + ha) Hall 2 


Dabei gilt für die zeitliche Ableitung H;(x) die Be- Durch Anwendung der genannten Beziehungen er- 


ziehung gibt sich nun im Wege partieller Integration: 
az a ie 
EC EZ H,(2) Bene (3). F,(2) =/f Jolt) » di 5 
R 
Die Funktionen H,(z) und H,(x) sind bei Wat- 2 : 
son tabelliert für reelle Argumente x zwischen Null 2:9 Fi). di 


und 15; sie befinden sich bi Jahnke-Emde J Pr 

im gleichen Intervall tabelliert, wo sie die Bezeichnung re dF,;(t) d 

8,(x) und S,(x) tragen. Die Ableitung der Gl.(2), die = Holt 

bei Watson in allgemeiner und daher etwas F 

weniger übersichtlicher Form gegeben ist, ist im An- 

hang kurz wiedergegeben. . =x.Fı(a) — j t+ Jo (t) di 
Ferner sei noch bemerkt, daß die Beziehung gilt: v 

© F,(&) = x Fı(®) — 2 Jı(@) 

FRLUL Wa T-, Seet (4). A = 

F,(x) = ft- F(t)-d — f t+ Jı(ty- di 


G.(4) folgt z.B. aus der bekannten Gleichung der 


z 
Laplace-Transformation = 2" P,(@)— SF, di+ 2 Jo(e) Filz) 


2-F,(&)=#'Fy(x) — Fı(z) + © Jo(®) 


L{Jy(at)}= a FIrtat)-e-P-dı 2-F,(2)=2 Fl) — Fu(2) —Fz(2)+ 2° Jı(®) 
? . (5) 3.F(a) = x Fs(&) — Fy(a) + © Jıl@) 


ar+i 
= 3-F =rc+»F u er 
Vote pH Fon en. 
naeın man setzt 4-Fs(2)= x Fila) — Fa(2) + Fılz) — x Jo(®) 
p=). j 4- Fe(z) = © Fl) — Fl) — Fu(z) + Fz(e) 
— 2: Jı(@) 
2. Mehrfach-Integration von Jy(z). 5-Fu(a)= © Fl) — File) + Fake) — 2 Jıle) 
Im nachfolgenden soll die Abkürzung verwendet 5-F,(2) = x Fale) — Frl) — Fs(2)+ Fs(e) 
werden: : — Fı(@) +2: Jo(e) 


6-F,(&) = 2° Fele)—Fy(&) + Fs(z)— Fı(2)+ 2° Jo(2) ‘ 


Fıla)=f Joldaı 
0 6-F,(&) = x + F,(&) — F;(x) — Fy(x) + F,(&) 


© x u 
Fa) = f Fr Wd=fdu-f Iolydi Er 
a a et 1, TFD= El) Fldt Rue) Fre)tee Ja) 
= 5 
F,(2) = F,()di = f IE .f Rd). di Das Ergebnis dieser Gleichungen läßt sich also 
fi ö v folgendermaßen zusammenfassen: 


a NE, Mehrfache Integrationen über die Bessel-Funkti 
= /S dvufduf Jylt)dt J,(x) lassen sich im ir der Rekursion ker 
070%, 120 führen auf Polynome enthaltend die Funktionen 


Also soll allgemein gelten: Jo(&), Jı(z) das Integral RAN) .di 
ö 


x 
Frzı (@) ec Pla (7).= sowie Potenzen des Argumentes x selber. 


3. Mehrfach-Integration von Jy(z) 
Für die weitere Rechnung werden die folgenden b- Geht man von der Gleichung 
kannten Gleichungen für Zylinder-Funktionen An- 


j dZy(x 
wendung finden: 2 ie = 2,-1(2) —Zp+ı(2) ». (11) 
5 i 
SeBWd=ahl) on... (8). @us, so folgt 
x E | 
an [Hd =f[2,0:d— 22,18) . (12). 
dz AusGl. (12) sieht ‚daßsi 
2 ( __ ren ea, (12)sieht man sich das ende Schema 


% 
| SA) dt = — Jul) 
Aus Gl. (9) ergibt sich im Wege ‘der partiellen ? h 


Integration: e 4 
; E ; LER =/Jl)d— 2. Jule) 
j ER Zott = 
[: Zi fa [HW-U- FAN 2: Ile) 
! : ’ 


7 
P (10). [ JA) di= IRAQ) dt — 2: Jy(e) 


= 2200) + [Zum 
% ot Man erkennt aus dem Schema auch, daß bei wei- 
ö teren Integrationen der auf den linken Seiten stehen- 


E = 
Kleine Mitteilungen 


den Integrale nur Auadenoke entstehen können, die > 


im Wege der Rekursion ebenfalls sich auf die Funk- z e ” 
 ionen” | = | fawal = \froal = [rad 
= en x SEE een 
Jo(&), Jı(z) und das Integral FELL, 0, 0 ne 6,0 0,7002 12,0 oma 
anne . y ’ ’ ’ Su —— , ’ 
. zurückführen lassen. 0,4 | 0,3947 | 6,4 | 0,7862 | 12,4 | 0,8106 
24, M ehrfach-Integration der Lom- a a a ee a 
i er “ E 7 x ’ ’ 3 ’ ’ j ’ 
ee ur ae is [ng ans 
: Führt man als Abkürzung für das erste er 1,4 1,1875 7,4 | 1,0720 | 13,4 | 1,0001 
_ über die Lommelsche Funktion 18 rn = a 2 an 
Au) = 2) 2,0 | 1,4258 | so | 1,2107 | 140 | 1.1907 
ı 2,2 | 1,4591 | 8,2 | 1,2402 | 14,2 | 1,1522 7° 
2,4 | 1,4702 | 8,4 | 1,2594 ‘| 14,4 | 1,1767 
> die Bezeichnung ein 2,6 1,4607 8,6 | 1,2679 | 14,6 | 1,1943_ 
Be 2,8 | 1,4323 | 8,8 | 1,2654 | 14,8 | 1,2041 
er ı(t 3,0 | 1,3876 | 9,0 | 1,2523 | 15,0 | 1,2052 
u 2 j I EM... (dd 30 | 19283 | 92 | 1.2292 | 16.0 | 1.1009 
2 3,4 | 1,2606 | 9,4 | 1,1980 | 17,0 | 0,9125 
3,6 | 1,1848 | 9,6 | 1,1593 |.ı8,0 | 0,8133 
2. entsprechend — um 3,8 | 1,1050 | 9,8 | 1,1153 | 19,0 | 0,8869 
4,0 | 1,0245 | 10,0 | 1,0670 | 20,0 | 1,0584 
E  Gntıla) = 7 nd... a) de | 09472 [102 | Loı2 | 21.0 | 1.1690 
“x 4,4 | 0,8751 | 10,4 | 0,9679 | 22,0 | 1,1224 
E20 = 4,6 | 0,8110 | 10,6 | 0,9207 | 23,0 | 0,9678 
er v, 4,8 | 0,7574 | 10,8 | 0,8774 | 24,0 | 0,8486 
Gay=2 | 0 = 5.0 | 0,7183 1:21,02) 0,8898 4. 2 
- 5,2 | 0,6864 | 11,2 | 0,8091 | © | 1,0000 
ee 5.42120,8716 111,4: | 0,7871 1 — 
& 5,6 | 0,6698 | 11,6 | 0,7734 | — — 
5,8 | 0,6819 | 11,8 — — 


0,7695. 


Ge) = | At): di 


0 
Bee ].i f* Ko} de 
nn | 
ya: = 7 = 
en (x) — 2 J Jı(t) di 
R or 


‚9,() +2: Sa) 
nn la) ie" -@ (2) — Gs(&) + 2° Fl) 
CHE CH er 
| =2Gs(e) — Pal nn?" "F en BET FEN 
E 2 “Fr, 2) 2 Bild1 
MER @ 2-F Zu 
rer ta * CE | Anhang: 
i Kae „leitung von G1.(2): Es gilt: 
ß toehilache Integrale Z le- Ju) H5@)] Sn 
nes BLOG EO oRe 
en le A) Hal ! a 
0020 hl) Hol) te Ji@)- Elta: hl- Biß) e 
Da = @)=—Jı 2)» ergibt sich durch Addition: 


eo). Aue) + A) Ha] = 


= * 2: Io): [a3 Dee Lu] + 
EN. CE o+ ne] 
J Re) : [eo 2) + z a 


= 2 "20 E 226) -1.% e]. 
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Durch Anwendung der Differential-Gleichungen: 
Be) + I Hife) + Hola) = — 
IE) + Io) + Jule) = 0 
folgt: » 
Eee): Hs) + 46) Hua)l) 


2 
= 2: Ju(e) F _ Ro) +2: Hol) ol) = + Jole) 


Also: f Jo) de= zZ [I o) He) + Jı(z) * Ho(z)] 


Berlin-Adlershof. G.B. Hagen. 


Zur Konstruktion der Achsen einer Ellipse. 
Das Problem, die Achsen einer Ellipse zu konstru- 
ieren, die durch ein Paar konjugierter Durchmesser ge- 
geben ist, wurde bereits von Pappus von Alexan- 
drien gelöst. Aus der Reihe der. im Laufe der Zeit an- 
gegebenen weiteren Konstruktionen, wird — seit un- 
gefähr 100 Jahren — die von dem Schweizer Mathe- 
matiker Rytz stammende Überlegung als ‚‚Stan- 
dartkonstruktion‘ verwendet. Im folgenden wird auf 
Grund rein kinematischer Betrachtungen eine andere 
einfache Konstruktion gezeigt und bewiesen. 


1.Konstruktion: MA und MC seien zwei 
konjugierte Ellipsen-Halbmesser (Bild 1). Man zeich- 
net durch C die Senkrechte 1 zu M A, trägt auf ihr die 
Strecke MA von C aus nach beiden Seiten ab und er- 


hält die Punkte EZ und F. Die Winkelhalbierenden . 


des Geradenpaares ME und MF sind die gesuchten 


Hauptachsen. Die Längen a und b der Halbachsen er- 
gebensich daraus, daß ME=a+bundMF =a— b 
ist. Das Halbieren des Winkels EMF läßt sich dadurch 
umgehen, daß man die Strecke MF auch auf dem 
anderen Winkelschenkel aufträgt: MF= MG. F@ 
ist dann parallel zureinen, also senkrecht zur anderen 
Achse. Außerdem wurde erreicht, daß GE = 2b ist. 
Der Halbierungspunkt H dieser Strecke liefert also 
die Halbachsenlängen MH =a, HE=b. 

2. Beweis: Werden die beiden als starr anzu- 
sehenden Strecken EK und FK (Bild 2) so bewegt, 
daß ihre Endpunkte Z und K (bzw. F und K) auf den 
festen Geraden ME und MB (bzw. MF und MB) 


gleiten, dann beschreibt der auf der Strecke EK (bzw. 
auf der verlängerten Strecke F.K) liegende Punkt C 


Bild 2 


die durch die gegebenen konjugierten Durchmesser 
AB und OD festgelegte Ellipse. Die beiden Polkurven 
für die Bewegung der Strecke EK sind der Kreis um 
M durch E und der Kreis über der Strecke ME. Der 
Radius ME der Rastpolbahn ist also gleich (a-+b). 
Das Analoge gilt für die Bewegung der Strecke 
FK-.MF=(a—b). Nun soll speziell jene Systemlage 
beider Strecken betrachtet werden, in denen sie jeweils 
auf den Geraden ME bzw. MF senkrecht stehen: 
E,K, bzw. F,K,. Aus den beiden kongruenten Drei- 
ecken (EK.M) und (K,Z,M) bzw. (FKM) und (K,F,M) 
folgt: MK=ME,=MF,. D.h. aber, daß die beiden 
betrachteten Systemlagen symmetrisch zur Halbieren- 
den des Winkels (EMF) sind. Der die Ellipse beschrei- 
bende Punkt 0 kommt dabei in die Lage C, bzw. Q. 
Dabei ist ZO=FO=EB,0,=F,(,. Die beiden Elli 
senpunkte C, und C sind somit ee 
trisch zur Winkelhalbierenden (EMF) gelegen. Diese 
Winkelhalbierende ist also eine Achse der erzeugten 
EHipse, , 
Braunschweig. 


R. Jakobi. 


BUCHBESPRECHUNGEN 


Dr. Wolfgang Gröbner (Prof. a. d. Universität Inns- 
bruck) und Dr. Nikolaus Hofreiter (Prof. a. d. Universi- 
en en Zweiter Teil. Bestimmte 

egrale. inger-V: i 
Innsbruck 1950. Preis Fe 

ie Vff. haben nun das im Kriege begonnene Unter- 
nehmen zum Abschluß gebracht a "grün lich 
durchgedachte und in allen Einzelzügen neu durch- 


gerechnete Tafel der Bestimmten Integrale heraus- 
gebracht. Der Band enthält solche Tote für die 
eine (unbestimmte) Integralfunktion nicht näher be- 
kannt ist, so daß nur bei spezieller Wahl der Grenzen 
Beziehungen zu anderweit bequem aufzubauenden Zahl- 
werten oder Funktionen eines Parameters bbar 
werden. ‚Der Aufbau der Tafel hieht nach einer 
Hierarchie der Integranden, die der im I. Teil über Un- 
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bestimmte Integrale entspricht. Durch Einführung ge- 
eigneter Parameter wurde versucht, möglichst viele 
Integrale gleicher Bauart zusammenfassend wiederzu- 
geben; doch wurde darauf geachtet, besonders wichtige 
Fälle ausdrücklich anzuführen. Die Verf. berichten über 
die Sicherungsmaßnahmen, die sie getroffen haben, um 
alle Angaben zu prüfen. Sie haben zudem, nach einem 
geeigneten System von Abkürzungen, zu jedem Integral 
einen Weg vermerkt, der erlaubt es zu berechnen; sie 
haben sich dabei auf eine Anzahl besonders wichtiger 
Methoden beschränkt, die eingangs kurz aufgeführt 
sind, nehmen aber keine Rücksicht auf Sonderwege, die 
vielleicht in einem Einzelfall rascher zum Ziele führen: 
das würde die Wegangaben zu sehr belasten. Die Wie- 
dergabe ist nach sorgfältig geschriebenen Formeln er- 
folgt; das Schriftbild ist vorzüglich. — Es sei gestattet 
zu vermerken: 1. In einigen Fällen wurden Abwei- 
chungen von den in der Literatur eingebürgerten Be- 
zeichnungen für nötig gehalten (Integral- und Diloga- 
rithmus, Eulersche Konstante u.a.); es sollte dann nicht 
bloß definiert, sondern auf die Art der Abweichung 
audrücklich hingewiesen werden. 2. Beim Integral- 
logarithmus wäre vielleicht ein Hinweis auf die grund- 
legende Rolle dieser Funktion bei doppelter Integration 
rationaler Funktionen am Platze (er fehlt auch in Bd.T). 
8. 8.174 Nr.421 finden wir mehrmals ‚‚Exponential- 
funktionen‘ an Stelle von ‚‚Allgemeiner Potenz“. So 
wird vom Zweig der Exponentialfunktionen gesprochen, 
der auf der positiv reellen Achse gleich dem Haupt- 
wert zu wählen sei; ebenso 8.178. 4. S.2, 011, 8 
wird das Symbol f(z)=g(x)-FO(x*) definiert durch 


lim Keen <C<o. Man kann das so machen; 


esist aber weder üblich noch auch zweckmäßig. Die 
Bachmann-Landausche Bezeichnungsweise mit den 
Großem O setzt nicht die Existenz eines Grenzwerts 
voraus. — Diese kleinen Erinnerungen werden aber 
der hingebenden und höchst gewissenhaften Arbeit der 
Verfasser keinen Abbruch tun, die wir gerne als vor- 
bildlich anerkennen. 


Gießen. Egon Ullrich. 


H. Grauden, Momenten-Einflußzahlen 
für Durchlaufträger mit beliebigen 
Stützweiten. IV + 90 Seiten mit 14 Abb. und 
80 Zahlentafeln. Berlin/Göttingen/Heidelberg 1951. 
Springer-Verlag. Preis 7,50 DM. 


Das neue Tabellenwerk zeichnet sich gegenüber den 
schon vorhandenen Tabellenbüchern dadurch aus, daß 
es auch für unregelmäßige Feldteilungen und wechselnde 
Trägheitsmomente schnell zu Ergebnissen von ausrei- 
chender Genauigkeit führt und dadurch überhaupt 
seine Berechtigung findet. Die Einflußzahlen der Ta- 
feln wurden mit Hilfe des Cross-Verfahrens berechnet. 


Dresden. Henn. 


Veröffentlichungen des Deutschen Aktuarvereins 
(Verein deutscher wissenschaftlicher und leitender 
praktischer Versicherungs- und Wirtschaftsmathe- 
matiker, e.V.) Nr.1, 538. Berlin-Charlottenburg 
1951. Robert Kiepert. 

Der nach dem Zusammenbruch für den Bereich 
Berlins wiedergegründete ‚‚Deutsche Aktuarverein“ 
nimmt mit dem vorliegenden Heft seine Veröffent- 


- lichungen wieder auf. Es enthält die folgenden Auf- 


sätze: B. Heinicke, Dr.phil. Eduard Rose zum 
Gedächtnis; K.Freudenberg, Betrachtungen 


über Sterbetafeln und Sterblichkeit; P. Lorenz, 


Der Schluß vom Teil aufs Ganze und der mathe- 
matisch-statistische Begriff der Wahrscheinlichkeit 
(Ein sinnentstellender Druckfehler findet sich auf 
Seite 27 Zeile 9 von unten. Es heißt ‚‚Formen ohne 
Inhalt‘ statt ‚‚Formen und Inhalt‘); J. Berger, 


= Über die Deckungsrückstellung der vertraglichen 


 Krankenversicherung. 


Dresden. Willers. 


Dr.R.v.Mises (Prof.a. d. Harvard-Universität Cam- 
bridge Mass. USA.) Wahrscheinlichkeit, Statistik 
und Wahrheit. Einführung in die neue Wahrschein- 
lichkeitslehre und ihre Anwendung. 3. neubearbeitete 
Aufl. IX + 2788. Wien 1951, Springer-Verlag. Preis 
brosch. 18,— DM. 

Das Wiedererscheinen dieses ausgezeichneten Buches, 
dessen erste und zweite Auflage in Bd.9 (1929) S.80/81 
bzw. Bd.16 (1936) 8.378 dieser Zeitschrift angezeigt 
wurde, wird allgemein mit Freude begrüßt werden. Die 
Grundgedanken der hier vertretenen Theorie, die beim 
ersten Erscheinen neu waren, sind heute mehr oder we- 
niger Gemeingut aller Autoren geworden, die sich mit 
den Grundlagen der Wahrscheinlichkeitsrechnung ke- 
fassen. Dadurch sind beträchtliche Teile der ersten Auf- 
lage, die polemische Auseinandersetzungen mit anderen 
Auffassungen enthielten, überflüssig geworden und 
konnten fortbleiben. Der dadurch besonders in den Vor- 
trägen 3—5 freigewordene Raum ist für Ergänzungen- 
verschiedener Art benutzt worden. So konnte Verf. auf 
die neueren Entwicklungen insbesondere in den angel- 
sächsischen Ländern eingehen. Dabei ergibt sich z.B. 
daß nach seiner Ansicht die von Fisher eingeführte 
likelihood und die daraus abgeleiteten Prüfverfahren, 
richtig angewandt, nicht über das Gebiet der hier ver- 
tretenen Theorie hinausführen, so daß Verf. die Hoff- 
nung aussprechen kann, daß dielikelihood bald wieder 
durch die einfachen und klaren Gedankengänge der 
Bayeschen Auffassung in dieihr zukommenden Grenzen 
verwiesen wird. Hervorgehoben zu werden verdient, 
daß die neueren Versuche zu einer Modifikation der 
Grundgedanken keine zwingenden Gesichtspunkte für 
eine wesentliche Veränderung der ursprünglichen Auf- 
fassung ergeben haben. 


Dresden. Willers. 


A.J. Chintschin, DreiPerlenderZahlen- 
theorie. (Herausgegeben vom Forschungsinstitut 
für Mathematik der Deutschen Akademie der Wissen- 
schaften zu Berlin.) 628. Berlin 1951. Akademie- 
Verlag. Preis geb. 6,50 DM. 

Der bekannte russische Forscher behandelt an drei 
Problemen, dem Satz von van der Waerden über 
arithmetische Progressionen, der Hypothese von 
Landau-Schnirelmann und dem Waringschen Problem, 
die modernen Methoden der additiven Zahlentheorie. 
Der elementare Charakter der Beweise, die keinerlei 
Voraussetzungen jenseits des Niveaus eines Abiturien- 
ten enthalten, wird dank der klaren Darstellung sehr 
gutherausgearbeitet. Daß elementar nicht identisch ist 
mit einfach — der Vf. betont es mehrfach —, zeigen 
alle drei vorgeführten Beweise. So schließt das Büch- 
lein mit der für die Zahlentheorie wohl in besonderem 
Maße zutreffenden Bemerkung: ‚‚Dieser durch seine 
elementaren Schlüsse wunderschöne Beweis wird 
Ihnen zweifellos sehr kompliziert erscheinen. Aber 
Sie brauchen an ihm nur 2 bis 3 Wochen zu arbeiten, 
um ihn vollständig zu verstehen und sich anzueignen. 
Gerade durch die Überwindung von Schwierigkeiten 
dieser Art wächst und entwickelt sich der Mathe- 
matiker.‘“ EHTS 

Als Einführung in die moderne additive Zahlen- 
theorie für weite Kreise von Liebhabern der Mathe- 
matik vortrefflich geeignet! Die deutsche Ausgabe 
muß darum dankbar begrüßt werden. 


Dresden. Draeger. 


Dr. Friedrich Hund (o. Professor a. d. Universität 
Jena), Einführungindie Theoretische 
Physik. BandI: Mechanik, 3. Aufl. 304 S. mit 
128 Abb. Band II: Elektrizität und Magnetismus. 
2. Aufl, 299 8. mit 75 Abb. Band III: Optik. 2.Aufl. 
2088. mit 100 Abb. Leipzig 1951. Bibliographisches 
Institut VEB. Preis geb. je 6,50DM. 

‘“ F. Hunds Einführung in .die Theoretische Physik 
ist in den letzten Bänden dieser Zeitschrift ausführlich 
besprochen worden, so daß eine neue Würdigung nicht 
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erforderlich ist. Das Buch hat seine alten Vorzüge be- 
halten, dieesin kürzester Zeitin die Hand aller unserer 
Studenten gebracht haben. Eine bemerkenswerte 
Änderung hat der Band ‚‚Elektrizität und Magnetis- 
mus‘ erfahren. Die Schreibweise der Gleichungen ist 
jetzt so gewählt, daß der Faktor 4 im Coulombschen 
Gesetz auftritt und aus den Maxwellschen Gleichungen 
verschwindet. Dies ist sehr zu begrüßen, da das Buch 
in dieser Hinsicht in unserer sonstigen Lehrbuch- 
literatur eine außenseitige Stellung hatte. 


Dresden. Recknagel. 


Dr. Kurt Reidemeister (o. Prof. a. d. Univ. Marburg), 
Einführung in die kombinatorische 
Topologie. (Die Wissenschaft herausgeg. von 
Prof. Dr. Wilhelm Westphal, Bd. 86), unveränder- 
ter Nachdruck. XII + 2098. Braunschweig 1951. 
Friedrich Vieweg u. Sohn. Preis geb. 12,— DM. 

Das Buch erfreute sich seit seinem Erscheinen 1932 
ständiger Beliebtheit und bedarf wohl kaum einer 
Empfehlung. Es wendet sich an Studenten der Ma- 
thematik der höheren Semester, die mit algebraischer 
Denkweise etwas vertraut sind, wiewohl es sachlich 
nichts voraussetzt. Dem Verfasser lagen vor allem 
die Zusammenhänge der Gruppentheorie und der 
Topologie am Herzen, da sich doch die Gruppen- 
theorie als das wichtigste Recheninstrument des To- 
pologen erwiesen hat. Deshalb sind die ersten drei 
Kapitel der Gruppentheorie gewidmet, insbesondere 
der Darstellung von Gruppen durch erzeugende und 
definierende Relationen und deren Bestimmung für 


Untergruppen. Die eigentlich topologischen Kapitel 
4—-6 können indes für sich gelesen werden, da auf die 
benötigten gruppentheoretischen Hilfsmittel jeweils 
sehr sorgfältig verwiesen ist. Das 7. Kapitel bringt 
verzweigte Überlagerungen in Anwendung von Grup- 
pen mit Operatoren. Esist etwas knapper gefaßt. 


Dresden Ott-Heinrich Keller. 


Colloque de Topologie (Espaces fibres), Tenu & 
Liege du 5 au 8 juin 11950, 137 S.,i Liege — |Paris 
1951. George Thone, Masson et Cie. Preis 175 francs 
belges, 1225 francs frangais. ! 

In einer gefaserten n-dimensionalen Mannigfaltig- 
keit geht durch jeden Punkt, einer r-dimensionale 
‚„‚Faser‘‘ (eine, Teilmannigfaltigkeit); alle diese Fasern 
sind homoiomorph; die Umgebung einer solchen Faser 
ist das ‚„‚kartesische Produkt‘ der Faser mit, einem 
(n—r)-dimensionalen Element. Die Betrachtung der 
gefaserten Mannigfaltigkeiten geht zurück auf die 
topologischen Untersuchungen von Poincare über die 
Richtungsfeldecr von Differentialgleichungen. Sie 
haben in den letzten Jahren ständig an Bedeutung 
gewonnen, vor allem durch ihre Anwendung auf Al- 
gebren und Liesche Gruppen. Der vorliegende Band 
ist der Bericht über eine Tagung vom 5.—8. Juni 1950 
in Lüttich, auf derH.Hopf,Eekmann (Zürich), 
Ehresmann,Koszul (Straßburg), Cartan, 
Leray (Paris, Hirsch (Gent) ihre Ergebnisse 
vorgetragen haben. 


Dresden, Ott-HeinrichKeller. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr.H.Schlichting (o. Prof. a. d. Techn. Hoch- 
schule Braunschweig), Grenzschicht-Theo- 
rie. (Wissenschaftliche Bücherei.) XV + 483 8, 
mit 295 Abb. und 32 Tabellen. Karlsruhe 1951. G. 
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